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Abstract. In this paper we derive a hierarchy of differential equations which uniquely determine the 
coefficients in the asymptotic expansion, for large A'^, of the logarithm of the partition function oi N X 
N Hermitian random matrices. These coefficients are generating functions for graphical enumeration on 
Riemann surfaces. The case that we particularly consider is for an underlying measure that differs from 
the Gaussian weight by a single monomial term of degree 2u. The coupling parameter for this term plays 
the role of the independent dynamical variable in the differential equations. From these equations one may 
deduce functional analytic characterizations of the coefficients in the asymptotic expansion. Moreover, this 
ode system can be solved recursively to explicitly construct these coefficients as functions of the coupling 
parameter. This analysis of the fine structure of the asymptotic coefficients can be extended to multiple 
coupling parameters and we present a limited illustration of this for the case of two parameters. 



The study of the Unitary Ensembles (UE) of random matrices I20j, begins with a family of probability 
measures on the space oi N x N Hermitian matrices. The measures are of the form 



where the function Vt is a scalar function, referred to as the potential of the external field, or simply the 
"external field" for short. Typically it is taken to be a polynomial, and written as follows: 



The partition function Zjy, which appears as a normalization factor in the UE measures, plays a central 
role in random matrix theory and its applications. It can be reduced to an integration over the eigenvalues 
which takes a form proportional to the integral below, for the particular case when k — N. 

When all the coefficients tk in the external field are set equal to zero the associated ensemble, corresponding 
to /iQ, is called the Gaussian Unitary Ensemble (GUE). Many simplifications occur in the Gaussian case (see 
[T^ for explanations of any unfamiliar terms): 

(1) The partition function, when all tk vanish, is a Gaussian integral, and can be evaluated exactly. 

(2) The matrix moments, / { Tr {M^)}''d^loiM), can be evaluated, using Wick's lemma, in terms of pair 
correlations of the matrix entries of M which are complex normal random variables. 

(3) The terms in these Wick coupling expansions are, in the manner of Fcynman diagrams, in 1-1 
correspondence with certain labelled, oriented graphs. 

These observations led to the conjecture |S1 that the logarithm of the partition function has an 
asymptotic expansion of the form: 



1. Motivation and Background 



dm = J- exp {-N Tr [Vt{M)]} dM, 
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where the coefficients eg(t) should be locally analytic functions of t. The Taylor coefficients of eg should 
enumerate topologically distinct labelled, connected oriented graphs that can be embedded into a Riemann 
surface of genus g in such a way that the complement of the graph in the surface is a disjoint union of 
contractible cells. Such a construction is referred to as a g-map (see section [131 for a precise definition). The 
eg{t) are generating functions for counting the number of (7-maps with given numbers of vertices of specified 
valence. This conjecture was proven in |2 for appropriate domains (see below). The present paper builds on 
these results to present a more detailed description of the coefficients eg(t) and related generating functions. 



More precisely, our interest is to develop a systematic, rigorous description of the fine structure for the 
large N asymptotics of the following family of integrals: 



(1.1) Z\^\ti,t2, . . . ,tv) 



exp { -N^ 



i=i j^e 

(1.2) V{X; ii,...,t,) = 14(A) = F(A) = iA2 + ^<,A^. 

where the parameters {ti, . . . ,ti,} are assumed to be such that the integral converges. For example, one may 
suppose that v is even, and > 0. We will sometimes refer to the following set of t = {ti, . . . , t^) for which 
(II. 1|) converges. For any given T > and 7 > 0, define 

v-l 

T(T,7) = {t e : |t| < T, > 7^ l^^l}- 

The parameter k is an integer that grows with N in such a way that lim^v^oo k/N = x, where a; is a finite 
non-zero value whose role will be specified more precisely later. 

In this paper we derive a hierarchy of differential equations which uniquely determine the coefficients in 
the asymptotic expansion of logZ^'' for monic even coupling parameters; i.e., we present the eg{t2v), for 
arbitrary v, as solutions to a system of ordinary differential equations. From this one can deduce functional 
analytic characterizations of these coefficients. Moreover, this ode system can be solved recursively in g to 
explicitly construct eg{t2u). We illustrate this process by constructing closed form expressions for eg{t2u), 
in which v appears as a parameter, for low values of g. This analysis of the fine structure of the Cg can be 
extended to multiple coupling parameters and we present a limited illustration of this for the case of two 
parameters: 63(^21/1 , t2i/2)- 

Remark In |2|, the so-called "Loop Equation" method is used to obtain some information about the fine 
structure of the coefficients. This approach is based on a formal derivation of a hierarchy of equations for the 
Cauchy transform of the mean density of eigenvalues. This interesting approach is unsatisfactory in that it 
relies on several interchanges of singular limits whose justification requires analytical considerations beyond 
the existence of the complete asymptotic expansion of the partition function. These analytical considerations 
are the subject of a forthcoming paper by Ercolani and McLaughlin ^21- D 

1.1. Leading Order Asymptotics. The leading order behavior of Z^\ti,t2, ■ ■ ■ ,tv) is rather classical, 
and is known for a very wide class of external fields V (see, for example, QHI)- We will require the following 
result. 

Theorem 1.1. There is Tq > and 70 > so that for all t G T(T,7), x e [1/2,1], and k/N ^ x as 
k, N ^ 00, the following holds true: 

(1) 

(1.3) Jim ^log{zP{ti,t2,...,t,)}^~I{x,ti,...,t,) 



where 

(1.4) I{x,ti, . . . ,i„) = inf 

Borel measures fi.fi>0,[ dfi—1 



- I vi\)d^,ix)- I I log\\~^,\d^ii\)d^l{r^) 



(2) There is a unique measure /iy which achieves the infimum defined on the right hand side of lil.4{ l- 
This measure is absolutely continuous with respect to Lebesgue measure, and 

duv — d\, 

= 2^X(o,/3)(A)V(A-o)(/3-A)/i(A), 

where h{X) is a polynomial of degree v — 2, which is strictly positive on the interval [a, /?] (recall that 
the external field V is a polynomial of degree v ). The polynomial h is defined by 



where the integral is taken on a circle containing (a, (3) and z in the interior, oriented counter- 
clockwise. 

(3) There exists a constant I, depending on V such that the following variational equations are satisfied 
by fiv : 

' 2\og\X-7]\-^dnv{v)+x^^VW > I for Xell\supp{^iv) 
(1.5) J 2\og\X^7]\-^d^iv{v)+x'^ViX) = I for Xe suppifiv). 

(4) The endpoints a and [3 are determined by the equations 

1-^ V'{s) 



a V(s - a){(i - s) 
^ sV'(s) 
a V(s - a){f3 - s) 



ds^O 
ds ~ 2'Kx. 



(5) The endpoints a(a;,t) and /3(x,t) are actually analytic functions oft and x, which possess smooth 
extensions to the closure of {x, t : x G [1/2, 1], t G T(T, 7)}. They also satisfy — a(l, 0) = /3(1, 0) = 2. 
In addition, the coefficients of the polynomial h{X) are also analytic functions oft andx, with smooth 
extensions to the closure of {x,t : x G [1/2, l],t G T(T, 7)}, with 

h{X,x = l,t = 0) = 1. 



Remark The variational problem appearing in (|1.4() is a fundamental component in the theory of random 
matrices, as well as integrable systems and approximation theory. It is well known, (see, for example, jMl)j 
that under general assumptions on V, the infimum is achieved at a unique measure fiy, called the equilibrium 
measure. For external fields V that are analytic in a neighborhood of the real axis, and with sufficient growth 
at 00, the equilibrium measure is supported on finitely many intervals, with density that is analytic on the 
interior of each interval, behaving at worst like a square root at each endpoint, (see ^U] and ^J). □ 

Remark We call the reader's attention to the parameter, x, in the formulation of the variational problem. 
We will consider the variational problem for x G (0, 1], and we are particularly interested in x near 1. This 
parameter represents the asymptotic ratio of fc to A^: x = limjv^oo k/N. □ 

Remark For a proof of H1.3(l , we refer the reader to ^H] > however this result is commonly known in the 
approximation theory literature. □ 

Remark It will prove useful to adapt the following alternative presentation for the function ij;: 
(1.6) ^(A) = -i,i?+(A)/i(A), AG (a,/3), 
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where the function i?(A) is defined viai?(A)^ = {X — a){X — P), with i?(A) analytic inC\[Q!,/3], and normalized 
so that i?(A) ~ A as A ^ oo. The subscript ± in R±{X) denotes the boundary value obtained from the upper 
(lower) half plane. □ 

1.2. Complete Asymptotic Expansion. In JS| it was established that a complete large N asymptotic 
expansion of II. II exists. In this paper we will use a straightforward generalization of this result: 

Theorem 1.2. There is T > and 7 > so that for t G T(T,7), and x — k/N in a neighborhood of x — 1, 
one has the N 00 asymptotic expansion 



The meaning of this expansion is: if you keep terms up to order k~ , the error term is bounded by Ck^ , 
where the constant C is independent of x and t for all t £ T(T,7) and for all x in the neighborhood of 1. 
For each j, the function ej{x,t) is an analytic function of the (complex) vector (x,t), in a neighborhood of 

(1,0). Moreover, the asymptotic expansion of derivatives o/log^Z^''^ may be calculated via term-by-term 
differentiation of the above series. □ 

Remark In 15 , this result was established in the case where a; = 1, and under the assumption that t G 
T(T, 7), for T small enough, and 7 large enough, so that Theorem 11.21 holds true. Under these assumptions. 
Theorem ^21 (below) was established. However, as observed in (Remark 2.1, page 2), the domain so 
defined is by no means the largest domain where the asymptotic expansion can be rigorously established. 
All that is required is the existence of a path through the space of parameters (values of x and t) connecting 
(a;,t) to (1,0) in such a way that all along the path, the associated equilibrium measure is supported on a 
single interval, with strict variational inequality on the support, strict positivity on the interval of support, 
and vanishing like a square root at both endpoints of the support. The collection of all such values of (x, t) 
defines a suitable candidate for a maximal domain, and the proof contained in 15 can easily be extended 
to show that the asymptotic expansion of the partition function holds on the interior of such a domain. In 
particular, the above Theorem may be easily deduced along these lines. □ 

Remark Recently, Bleher and Its 8 have carried out a similar asymptotic expansion of the partition 
function for a 1-parameter family of external fields. A very interesting aspect of their work is that they 
establish the nature of the asymptotic expansion of the partition through a critical phase transition. □ 

1.3. Graphical Enumeration and the Partition Function Expansion. Our goal in the work we present 
here is to establish analytical characterizations of the coefficients Cg and, when possible, to derive explicit 
expressions for these coefficients. This is what we mean by the fine structure of the expansion. 

In addition to providing the first proof of the asymptotic expansion described in Theorem 11.21 [TH] also 
provides a very detailed explanation of the connection between the asymptotic expansion and enumerative 
geometry, originally investigated by physicists in the 70s and 80s (see, for example, [5], JSli and references 
contained therein). Equipped with the existence of the asymptotic expansion (and the subsequent result 
that it may be differentiated term by term) , one shows that there is a geometric characterization of each eg 
as a generating function for enumerating topologically distinct embeddings of graphs into Riemann surfaces 
of genus g. 

A map D on a. compact, oriented connected surface X is a pair D = {K{D), [^]) where 

(1) K{D) is a connected 1-complex; 

(2) [i\ is an isotopical class of inclusions i : K{D) — > X\ 

(3) the complement of K{D) in A is a disjoint union of open cells (faces); 

(4) the complement of Ko{D) (vertices) in K{D) is a disjoint union of open segments (edges). 

The Cg enumerate labelled maps. To be precise we introduce the notion of a g-map which is a map in which 
the surface X is the closed, oriented Riemann surface of genus g and which in addition carries a labelling 
(ordering) of the vertices. 



(1.7) 
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Theorem 1.3. |15| The coefficients in the asymptotic expansion j i. ?| ) satisfy the following relations. Let g 
be a nonnegative integer. Then 

(1.8) eg(ti...t,)= V — i -(-ti)"^ ...(-t,)"-Kg(ni,...,n,) 

^ nil . ..nj 

nj>l 

in which each of the coefficients Kg{ni, . . . ,n„) is the number of g-maps with Uj j-valent vertices for j = 
l,...,v. 

1.4. Outline. The organization of tliis paper is as follows: In section 2 we present the new results concerning 
the fine structure of the Cg and related generating functions that will be proven and further explained in the 
remainder of the paper. 

Section 3 is concerned with the leading order term, cq. The results here are fundamental for the charac- 
terization of all the higher order terms. We derive closed form expressions for eg as a function of each of the 
valence coupling parameters t2u . We also relate these evaluations directly and explicitly to the enumeration 
of planar graphs. 

In section 4 a continuum limit of the Toda Lattice hierarchy is rigorously derived in which the hierarchy 
of Toda times corresponds to the valence coupling parameters t2u- This continuum limit is then used to 
derive another hierarchy of differential equations whose solutions are the Cg. 

Finally in section 5 we show how the differential equations derived in the previous section are used to 
inductively generate explicit expressions for the Cg. From this we characterize the function-theoretic structure 
of the eg as well as present explicit formulae for the Cg for low values of g. We also show how our results 
may be extended to the case of multiple times. 



2. Results 

For eo we have explicit formulas for monic even times 

Theorem 2.1. For potentials V of the form V ~ + t2iyX^'^ , the asymptotic expansion p.7|l holds true 
for all t2u > 0, and in addition, we have the explicit formula 

eo = T]{z- 1)(2- r) + ilog(z) 

where 



4 ■ 



Here 4z can be interpreted as the global analytic continuation of 0^ which determines the support (— /3, 13) of 
the equilibrium measure. The variable z is locally an analytic function of t2u, which satisfies the algebraic 
relation 

1^ z + 2iy{ _^ \x''-^t2uz'' . 
The singularities of cq occur at z = and z = oo. The time derivative 

is polynomial in z. 



One also has a local analytical representation (here the index n2v is replaced by n, so that ko{Q,. . 0, n2v) 
becomes Ko{n) ), 



Ko{n) = {cu) 
Cu = 2v 



(i/n- 1)! 



((z/-l)n + 2)!' 

2v- 1 
v-1 

where KqIu) ~ HQ{n2,y) is the generating function for 2v-valent 0-maps. 

To get a handle on how the higher coefficients Cg depend on the parameters t = t2v we exploit a remarkable 



relation between the partition function zjj^^(t) and the solutions to the hierarchy of completely integrable 



semi-infinite Toda lattice equations. These differential equations may be succinctly expressed through the 
semi-infinite tri-diagonal matrix 



(2.1) 



/ 





1 





... \ 




bl 


1 
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bl 





The Toda Lattice system at level 2v can then be defined as 



(2.2) 
(2.3) 



2 d£, 



k+2,k 



fc+l,/c-l J 



(£^'')fe+i,fe-i = 

ji,»2,---,«2i.+i;Nj+i-»j l=i;»i='=+i'»2i-+i=fe-i 
The sum here is indexed by walks of length 2v along the ID integer lattice from fc -|- 1 to A; — 1. The solution 
of this system may be expressed directly in terms of the partition function Zk{t\,t) = z'^^ associated to the 
potential V = ^A^ -|- iiA -h tX'^": 



(2.4) 



1 

bliO = fc(^)^log^fc(ii,s)t^^_fe-i/25^=o,5=2Cfe- 



As a dynamical system, one is really considering an initial value problem, with 
(2.5) bkiOf = k. 

We can now state our next main result which characterizes the continuum limit of the Toda lattice hierarchy. 
Theorem 2.2. For all t>Q,h\ has a valid asymptotic expansion of the form 



1 



1 



&fc ~ k{zo{s) + -^ziis) + -^Z2{s) + ■■■) 



A:2 



where s = —2k" t. The terms of this expansion are determined by the following partial differential scheme: 



where 

fis,w) : 



fo{s,w) + -^fi{s,w) +■ 



evaluated at w=l 



f{s, 1) = zo{s) + ^zi{s) + ^Z2{s) + 

fg{S,W)=W^-^Hg{w''-^S). 



Note that bf, and kf(s, 1) possess the same asymptotic expansion. 

The forcing term F^^\- ■ ■ )\w=i is a homogeneous multi-nomial of degree v + 1 in the f^M which does not 
contain any instances of Za for a > j . 

These forcing terms have the following form: 



2g+l , ^ r,{V) 

V:\V\=2g+13 p(V)<iy+l j=l 

where V = \J'JZ\ Vm is a partition of2g + 1; rj is the number of times a "part", Vm, of cardinality \Vm\ = j 
appears in the partition; p — p{V) = '^f'j{V); and 



d 



■ ' <2i/ 



Llj=l ' J- i<ii<.<ip( 



i^coeff ofx'^-P'^^^+^in P{x) - coeff of x'^'p'^^^^'^ m Q(x)) ; where, 



p{V) / m PI \ 



V ' dxs 

a^S m — l \ s — 1 



(1 + xif'-' ■■■{! + x,r=~'-i-i • • • (1 + xpy^-'^-'-\i + xp+i)2-v 1^^^^, 

d 

dxs 



p(y) / m f) \ 



a^S m—l \ s—1 

• (1 + xir^-i • • • (1 + x,r=-'»--i • • • (1 + a;p+i)'''-'-u,=.. 

Remark We refer to the above as a partial differential scheme because it signifies not an equation to 
be solved but rather a prescription for generating a hierarchy of ordinary differential equations for the Zg. 
The ode hierarchy is constructed from the scheme as follows. The ode at level g is obtained by replacing 
the expansion / and its w-derivatives, fwifww, etc., in the pde scheme by their order g truncations. Then, 
equating the coefficients of all terms of order k~^^ in this truncated scheme and setting w = \ yields a 
fc- independent ode in s which is the g^^ equation of the continuum- Toda hierarchy. This is an ode for Zg{s) 
in terms of Zj{s) for j < g. □ 

The first equation in the hierarchy (i.e., the one coming from the fe'^-coefficients of the above scheme) is 
a nonlinear ODE for zq: 

with the initial condition zo(0) = 1. This ODE is solved implicitly by a solution to the algebraic equation 

(2.6) l-zo(.s)-c,szo(s)". 

As indicated in Theorem ll.il (5). this relation can also be derived directly from the characterization of the 
equilibrium measure. Relation H2.6|l allows us to write s as well as derivatives of zq as rational functions of 
zq. This will be exploited to arrive at the explicit representations of the Zg given below. 
The equation in the hierarchy is linear in Zg, and can be written as: 

(2.7) z'g{s) = {fl^fg^ + i^fr'faUu.)^^, + ForcingJ^.i 
where 



(2.8) Forcing = 



<i-j < g 

\ ii-\ h i^+i = g J 



F^-^ [2g - 2] + F^"^ [2,9 - 4] + ... + F^^) [0] , 



and Fj^'^\2r] denotes the coefiicient of k in F^^K We note that the terms in Forcinggj^^i depend only on 
Zj,j < g and their derivatives. 



Remark Amongst the results on the above singular limit of the Toda lattice in the literature, we remark 
that the recent work of Bloch, Uribe, and Golse IB) is related, in the sense that in their work, through the 
use of the theory of Toeplitz operators, the authors establish the existence of an asymptotic expansion for a 
continuum limit of a finite dimensional Toda lattice. □ 

For the case of planar maps {g = 0) such generating functions have received significant attention recently 
[7||^. Our explicit calculations presented later can provide a basis for extending these studies. 

Before proceeding to the statement of our next result we need to introduce a scale of function classes that 
will enable us to describe the functional nature of the coefficients Zg that we have just introduced as well as 
that of the generating functions eg. We will refer to the classes as iterated integrals of rational functions or 
iir for short. These classes are defined inductively in terms of the variable z = zq regarded as an independent 
variable as follows. To begin with, the class contains rational functions of z. One then adds integrals of 
these rational functions with respect to dz. Next one considers the vector space of polynomials in products 
of these integrals over the field of rational functions in z and augments the space by integrals, with respect to 
dz of these functions. Then take the vector space of polynomials in these latter integrals and add integrals of 
these. One continues this iterative process up to any given finite stage. These are the classes of functions we 
refer to as iir. In our case the rational functions at any stage will be restricted to the sub-ring of functions 
whose poles are located at either z = or 2: = i^/{>^ — 1) for a fixed value of v. These classes of functions 
will certainly include the class generated by polylogarithms 17 but may be larger. 

Theorem 2.3. (1) The coefficient Zg is of class iir in zq with singularities only possible at zq — and 
zq = v/{v -I). 

(2) The coefficient Zg is more explicitly presented as a function of zq through the following integral 
solution of Equation {2. ?| ); 

Zo{sf^^-<^^ p(-) {v-{v -l)y)^ 
V - {v - l)zo{s) J-^ Cuy" 



= .. / „ .,.+3-2g Forcingg(y)dj/. 



(3) In the above equation, Forcing^, formerly a function of a great many arguments, is in fact a function 
of zq{s) alone, which will henceforth be denoted as Forcingg(zo). 

We also derive a hierarchy of differential equations for the eg{t2u) with data given in terms of the z'^s. 

Theorem 2.4. The g'th equation in the hierarchy of equations governing eg{t2p) is 



(2-9) 1^ K-'^e. (-''"^^)] 



9 9 /^(2n+2) r 

V = - u;2-2(ff-»)e i-w^'-^s^ 



(271 + 2)! aw(2"+2) 
the fc^^f term of log ( X! ) • 

\n=0 / 



Equation H2.9|l determines eg(— s) from a second order differential equation for Cg with forcing terms 
depending on e„, n < g, Zn, n < g, and their derivatives. 



Remark Observe that the RHS of H2.9|) is a function of zo{s) which will henceforth be denoted by 
driverSg(zo)- D 

Theorem 2.5. (1) The coefficient eg(— s) is of class iir in zq with singularities restricted to zq = and 
Zq :=^v/{v -I). 
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(2) The solution of \2.y\) may he represented as 



(2.10) e,(-.) = -^-__i-__^drivers,(zoW) 



2-25 V^o(s) - 1/ Ji \ CuU 

(l-2g) Uo(s)-lj ii U-y 

;f^s(29-2)/(i^-l) ^ ^^g(2ff-l)/(^-l) 



(driversg(y))*dy 
(driversg(2/))' 



when g 1, where Ki and K2 are constants of integration either determined by the requirement 
that Cg be a locally analytic function of s or by the evaluation of Cg for low values of v through its 
combinatorial characterization; and, when 5 = 1, 



ei(-s) = 



driversi(y)(iy 



(2.11) =-^logiiy-{iy-l)zo{s)), 

where we have chosen the principal branch of the logarithm. By (driverSg(?/))* we mean the derivative 
0/ driverSg(2/) with respect to y. 

3. Leading Order 

We will show that the leading order coefficient, ea{x, t), of the asymptotic partition function, fc^^ log(Z^-' {t)/z'^^^ (0)) 
is found in terms of the equilibrium measure /i = IJ-Vt/x- Define i^iX) by 

dfi = 'tjj{X)dX. 

We denote the leading order behavior 



Et = lim — 



^ log [Z^^\t)) =I{x,t), 



and 



(3.1) lim -7^ log -^Y^ l^Et- Eq. 



1 ..J zPit) 



More explicitly, by (|1.4|l 



Et = / ]Md^(A)+ / / log|A-r7r'rfM(A)dA^(r/) 



X 



(3.2) = (^-^,V(A)) -(£V», 
where 

(£/) (A)= I log|A-ry|/(77)dry 
is the logarithmic potential of the measure f{r])dri, and where the inner product (•, •) is defined by 

if, 9)^ I f{X)9iX)dX. 
Using H3.1|l and \d.2\i together with (|1.7|l . we find that 

(3.3) eo(i, 2;) = + ^A2^ ^(A)") + (r^, V^) + i Qa^, V.o(A)') - (£^0, ^0) • 



We recall here the parameter I mtroduced in p.5|l : 

21og|A-s|-V(s)'^s + ^^^ 



= -2(/:7/^,v^) + ^. 

X 

Using this we have the following reduced formula for eg: 

(3.4) eo = + _ + _(A , V^o) - 

This formula shows that there are two fundamental quantities that need to be calculated in order to evaluate 
cq. These are the moment {V,'>p) and the lagrange multiplier I associated to the constraint that the measure 
fjL should have total mass — 1. The other quantities appearing in (13.41) . (A^,?/'o) and Zo are evaluated by 
specializing the fundamental quantities at i = 0. Evaluating the fundamental quantities will require an 
explicit asymptotic expansion, for large A, of the equilibrium measure 'ip{X)d\. We develop this in the next 
section. 

3.1. Explicit asymptotic expansion of the equilibrium measure. The equilibrium measure fi is of 
the form 

(3.5) V(A) = ^h{X)y/>?^ x[-mW. 
where /i is a polynomial determined by 

(3.6) MA) = , ^ "^'^"l -ds, 

the integral in 1)3. 6|l is taken along a simple closed counterclockwise contour large enough to contain the 
interval [— /3, [3] and the point A. 

We will now evaluate the polynomial h by expanding the integrand of (|3.6|l for large A and calculating 
the loop integral on this expansion. 

Define the sequence {vj}J^Q by 

1=0 

whose Taylor coefficients can be computed to be 

with defined to be 0^/2. 

Expand the polynomial /i, given by (|3.t)|) . in terms of its coefficients hj: 



(3.9) h{\) 



Next note that 

2{i + l)v^ 1 



(3-10) =^ ■ ,n 

A direct computation of (|3.6|l gives 

(3.11) M^4M^-j>^^. 
The constraint that the total mass of /i should be 1, 



(3.12) 1 = / V(A)dA = -(p V(A)dA, 



f 1 
^{\)d\ = 



10 



implicitly determines i as a function of 0^ . Here, the integral is over a contour containing the interval [— /3, /?]. 

We compute the loop integral in (|3.12() over a large contour and find that the constraint can be expressed 
as a relation between the coefficients {/ij} and {vj\: 



(3.13) 2 = !^+gM. 



X '■ — ' X 

J=0 



The relation H3.13|l simplifies to the expression (using identities (13 .81) and (|3.11|l 



If we set z = /3^/(4x) H3.14|l becomes 

(3.15) l = z + 2i^|^"' ^y-Hz"". 
When 1/ = 2, 

(3.16) /i(A) = i(l + 2t/32+4tA^), 

X 

with the constraint 

(3.17) \ = z^Ylxtz'^. 

3.2. Explicit calculation of I. We will first derive an analytic expression for I in terms of /3^. To this end 
we study the logarithmic potential of [i: 

(3.18) g{\) = / \og{\~ s)i;{s)ds. 

J-p 

This function is analytic in C \ (— oo,/3]. For A G (— cxd,/3] we define two functions (7+ and g- by 

g±(A) = linig(A±ie). 

Choosing the appropriate branch of the logarithm we find that these functions are expressible as 

(3.19) 5±(A) = j \og\\ - s\i!{s)ds±i'K j'^ ^{s)ds. 
A calculation using H3.19|l and (|1.5|l shows that 

(3.20) .g+(A) +.g_(A) - + / = 2 / log |A - s|V(s)ds - + / = 0, 

X J x 

and 

(3.21) 5+(A)-g_(A) = -2z^ / V^(s)ds. 
Equation H3.21|) implies that 

therefore, if A is in C/(— oo,/?], then equation (|3.2()|l becomes 

2g{X)-2iTr [ i;{s)ds - + I = 0, 

and 

(3.22) -l^2g{X}-^^^+ [ h{s) ^/ - (i^ ds . 



This is the basic expression we shall use to calculate I. Since / is a constant it can be evaluated for any 
choice of A. We will evaluate it by studying the limit of (|3.22() as A —> oo. First observe that in this limit 

(3.23) 5(A)=log(A) + oQ 
Thus the principal issue is to develop an asymptotic expansion of the indefinite integral 

/•A 

(3.24) / h{s)^/s^ - f3^ds. 
The details of the derivation are deferred to Appendix^ the result is the following: 

(3.25) |\(.)V^^ = iwiX){X'-p')^^^^ + ^X'^^^p^ 



-2 log 



A , v/A^ 



■2 



[3 (3 



where 



(3.26) W(A) = ^WpA2p, 

p=i 



with 



1 "'^ 

(3.27) Wr, ~ - — — > . 

^ ' ^ 2vJp+l) ^ ^ ^ 



3=P 



As a result of these calculations and taking the limit as A ^ oo, we deduce the formula 



(3.28) -I = w,v,.^ ~ J2 wjv, - 1 + log 

in which terms of non-constant order (which must cancel in any case) have been dropped. 

Remark The relations generated by setting the non-constant terms equal to zero are equivalent to the 
moment conditions ll.il f4'l for the measure (see also jH]). □ 
Equation H3.28(l simplifies to the expression: 



(3.29) _,^4(2. 1)(2. 3)!.g:_ 

^ ' xv{v- 2)!2 4^ V 4 



Next, we use (j3.14|l to express t as a function of 0^ and 



iv - lV24''-i 

("») '-k^'"--^'' 



Substituting (|3.3U|I into (|3.29|) we finally have: 

As an example, consider v = 2 for which equation (|3.31|l becomes 



3.3. Explicit calculation of (VjV')- Next we need to find an expression for {V,ip): 



(3.32) 



Evaluating l|3.32(l by regarding its double as a contour integral and computing the value in terms of Vj and 
hj, we find that 

(3.33) (V, ^) = ^ g ^^.-.+1 + £ + i E + 

j=0 3=0 

This expression simplifies to 

4-1-^/3^-+^ (2^)! 24-^-^^(2^ ~l)P/34- 4-i-(^ + l)t/?^'-+^(2z.)! 

^ ^ + !)(;.- 1)!2 ^ 32a; - 1)!4 ^ a;(i. + l)!2 

As before, we simplify equation (|3.34l) by substituting the expression H3.30|l for t in terms of /S^ : 

3.35) V') = 5:77 — — TT • 

When = 2, we find that 

3.4. The explicit formula for eo{t). We can now put all of the pieces together. Evaluating the above 
expressions when t = (or equivalently when /3 — 2^/x), we find Eq — — 3/4 + 1/2 log(x). Collecting all the 



components ] 


gives 






(3.37) 




+ ilog( 


Ax) 


where 








(3.38) 


V 




1)2 




+ 1)^ 


(3.39) 


r 


3(j^4 
v — 


-1) 
1 ' 


and 








(3.40) 




4x 





We have proven that cq has an explicit representation depending only on 1/ (see the first part of Theorem 
lOl . 

When 1^ — 2 the result is 
(3.41) eo = ^(z-l)(z-9) + ilog(z). 



Remark It is interesting to compare our expression here to a formula derived in 8 . The difference is 
that while our formula is more explicit, it is restricted to a 1-parameter family of times corresponding to a 
fixed valence of the vertices. □ 
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3.5. Enumeration of Planar Graphs. It follows from Theorem 1 1 . 31 that the coefficients, (i) , of the 
leading order term eo{t) count the number of planar 2j/-regular maps with j vertices. We have shown that 
eo(t) can be explicitly expressed in terms of the auxiliary function z{t). The latter solves the algebraic 
relation H3.15|l . So to find the Taylor Coefficients of eo{t), we must first find the Taylor Coefficients of z{t). 
We define 

(3.42) a = -c^x^-h, 
where 

Then the polynomial relation defining z(t) is 

(3.43) l^z{t)-az{tY 

in a neighborhood of t = 0. (We view x here as a scaling parameter. When x is set to 1, we will recover the 
counting function.) 

The j*'' coefficient of the Taylor expansion of z as a function of a near 0, z = J2j>o Q'^'' ^ '^^ course 
given by 



1 / -2 a , 



Making the substitution u — z(a), and using the evaluation 

dz z" 

(3.44) 



da 1 — vaz'^ ^ 

derived by differentiating the relation H3.43|l . and also using H3.43|l to eliminate a, this integral becomes 

1 / (v- (v- l)u) u"^ , 
3.45 — i ^ du. 

Applying the binomial expansion, it is then straightforward to evaluate this loop integral and find that 

We note that these coefficients are precisely the higher Catalan numbers, which play a role in a wide 
variety of combinatorial problems. For a discussion of these applications and their relation to the work 
discussed here see '22^ . 

A similar approach will yield coefficients for the other terms in the expression H3.37|l for eg. For instance 
the Taylor coefficients of log(z(a)) = 12'jLi ^j"''^"' given by 

_(^) _ 1 / log(M)(l - i/au''-i)uJ'''_ 



Expanding the integrand in the vicinity of u = 1, we can evaluate these coefficients as 
(3.47) Lf^= ^ ' "'-^^^^'^^^ - ^ ^ "^-^ 



By the same method one also derives the coefficients in the expansion of the quadratic term, {z{a) — 1)^ 



(3.48) 



(^) _ 2 / lyj 
J \J - 2 
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Finally one may assemble all these contributions in (|3.37|) and substitute for a as in 13.42|l to conclude that 
the Taylor coefficients of cq with respect to ^21/ satisfy 



-{r - 1)7] 



2(j-l)(i.j)! ^ 1 



After simplifying this equation we get 
(3.49) .0 
This establishes the second half of Theorem 12. II 



" ^^'^ "{{,^-111+ 2)1- 



4. Continuum Toda Equations 

We will show that the Toda Lattice equations at level 2j/, with the initial conditions bk = Vk, possess a 
continuum limit under an appropriate scaling. These equations will lead to a description of the evolution 
of the asymptotic partition function coefficients, &g{t). More precisely, in Section 4.1 we establish Theorem 
2.2 and in subsection 4.2 we establish Theorem 2.7. To this end, we first review some basic facts about the 
Toda lattice equations. 

The Toda Lattice Hierarchy is given by the differential equations 

(4.1) ^ = [i3.(L),L], 

where L is a symmetric tridiagonal matrix of the form 



(4.2) 



L = 



/ ai 


bi 











...\ 


bi 




b2 













b2 


03 


^^3 













^^3 




bi 














/ 



and Bi{L) ~ (^')+ ~ (-^*)- '^h^re the plus subscript denotes upper triangular projection and the minus 
subscript denotes lower triangular projection. 

We briefly recall the notion of tau Junctions in terms of which the solution of the Toda Lattice equations 
at different time parameters can be expressed. Let Lq be a semi-infinite tri-diagonal matrix of the type that 
arises as a 3-term recurrence relation for orthogonal polynomials with exponential weight. The dependence 
of the tau function on the is given by 

(4.3) - '^ ^ • ^-^ - '-^"^ rA.r2ni/2 



,(6,...,^.;io) = (det [M]]) 



where M — ey^Y>(Yli"i=i ^i^h) Mj is the j x j upper left block of M in the basis of orthogonal polynomials 
with respect to d/j,. From the matrix factorization method of solving the Toda equations ^H! one finds that 



(4.4) 



Furthermore, from 14.4|l and the differential equations for the ^i-Toda flow, one can deduce that 
(4.5) 



The representations (|4.4() and H4.5() of bj in terms of tau functions are referred to as Hirota relations, familiar 
from the theory of soliton equations. 
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Finally, through the use of Hankel determinants, the tau function may be related to the partition function, 
defined as follows 



Zl^' . Recall that the Hankel determinant associated to the weight dLu{X) = exp(^^^j^ ^iA*)dA is 



(4.6) 



Co 

Cl 



Cl 
C2 



Cfe-1 
Ck 



Cfe-1 Cfe • • • C2fe_2 

where Cj = J^X^ duj{X) are the moments of dio. There follows the relation 

The partition function may be obtained from the Hankel determinant by a rescaling of the time parameters 
(see, for example, jEl Equation (1.29)]). Consequently one can relate the partition function to the tau 
function, in the case of interest to us, as 



(4.7) 



ZkiO) 

where Lq is the semi- infinite tridiagonal matrix with 6fe — ^/k, = 0, and where ^ = S^2u- For this initial 
condition, the evolution of L under the even 2i/-flow preserves the vanishing of Ofe. It will be natural to study 
this Toda evolution in terms of the self-similar variable 

(4.8) s = 2k''-^^. 
We will start from the wnscaZed partition function ()l.l|l : 

(4.9) zf)(_^,,_^) = 



/.../exp|-EQ^-?^"-^i 



A, } V(A) d'^X, 



where V(A) = ]^|Aj - A^]. 

Using the Hirota relations, 1)4. 4|l and 1)4. 5|l . we have the following Proposition. 
Proposition 4.1. The Toda recurrence coejficients may be represented as follows. 

1 d2 



(4.10) 



log [Zk{h,-s)] 



where s < 0. 



ti=0 



2k dti"^ 

In addition, the following asymptotic expansion holds true. 

(4-11) ^I = ^E^^'W^' 

where Zg{s) is an analytic function of s in a neighborhood of s ^ 0, and 

(4.12) zo(0) = 1, ^,,(0) = 0, 5 > 0. 

To prove this proposition we scale the eigenvalues in the unsealed partition function as Xj = N^/'^Xj under 
which 



^(fe) ^ ^feV2 



/^../o,p(-»{|:(iA3-«A.-sr--«LA,)}] 



V(A) rf'^A 



(4.13) 



^fe-/2^(fe) _ ^^^.-1 
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We observe that the Toda variables associated to our random matrix ensemble with the corresponding N 
scaling have the following representation: 



(4.14) 



bl 



2dCi 
2N dti 



^ log Zl 



2 log 



(k) 
(k) 



■'N 



N' 



■log 



-e,N 



1~ % 



We next define a spatial scale through the ratio between k and N: set k = xN . If wc furthermore 
introduce the scaling \j = \/xXj , then we can rewrite 



(4.15) 



/ . . . / exp |- Je (^A,^ + i--'^>^r + ^A,) I j V(A) d'^X 



k^/2y{k) ( h 



Putting all the above together, the natural variable to consider for this continuum limit of the Toda 
Lattice Equations is 



(4.16) 



1 (k) 

One may conclude directly from the main theorem of 15 that p- log Zf. has a full asymptotic expansion 
in powers of uniformly valid in an appropriate {ti, — s) domain which includes arbitrary negative s and 
simultaneously an open neighborhood of =0. On this domain the asymptotic expansion of derivatives 
may be calculated via term by term differentiation of the original asymptotic expansion. Moreover, the 
coefRcients in this expansion are complex analytic in a neighborhood of (0, 0). It follows that 



(4.17) 



1 d^ 
2fc2 dh^ 



log (ti,-s) 



has a full asymptotic expansion in powers of fc ^ uniformly valid for negative s and whose coefficients are 
analytic in a complex neighborhood of s = 0. 



4.1. The Continuum Toda Hierarchy. In this subsection we prove Theorem 2.2. 

We begin with the Toda lattice systems (14. If) introduced earlier. It will prove convenient to transform 
this to the so-called Kostant-Toda lattice form 14| where L is asymmetric tri-diagonal but in which the &^ 
are the fundamental variables. This is achieved by conjugating L in H4.1|) by the diagonal matrix 



(4.18) 



D = 



/I 












^ib2) 










The fundamental matrix variable for the Kostant-Toda system is 



(4.19) 





1 
















1 
















1 











bl 





1 



V 



Because D is diagonal it is straightforward to see that the RHS of H4.1|) transforms as 

(4.20) D-^[B24L),L]D = [D-^B2^{L)D,D-^LD] = [B2^{C),£]. 
One also has 

= (D-^LD)^ = D'^L^D +[C,D~^D^]. 
Combining this observation with (|4.2U|) in H4.1|l we see that the Lax equation transforms to 

(4.21) C^ + [D-'D^,C]^[B2A^),C]. 
By direct calculation one can check that 

where 



(4.22) 



Using this observation the Lax equation H4.21|l becomes 

1 













/O 1 











...\ 





1 
















1 
















1 




V : 








/ 



(4.23) 



--C^^[B2An,J^]-[D-'D^,e] 



Since the LHS is lower triangular, we may reduce this equation to 

1 



(4.24) 



-£5 = [S2.(/:),£] 



Since C is tri-diagonal, the RHS may be further simplified as 

[i?2.(/:), C]_ = {C?^ - £2-, £] _ = - {C?^, £] _ . 

Finally, since L = £_ + e with £_ supported only on the first lower diagonal, it follows that £^ is only 
supported on the first lower diagonal and therefore we may use the following final form of our Lax equation: 



(4.25) 



1 



2-'. L---J(-i) 

where the subscript (—1) denotes projection onto the first lower subdiagonal. At the level of the matrix 
entries h\ this system of equations becomes 

(4.26) 

where, 
(4.27) 



(£^'^)fe+i,fe-i 



)fc+2,fc — (£ 



ii,i2,---,«2(y + i;|ij + l— ij I — A:+l,i2iv + i— fc — 1 

as stated in (|2.2|l and (|2.3|) . 

Note that in (|4.27|) the sum may be viewed as being taken over walks, w, from k+1 to k-1 of length 2^. We 
will sometimes represent a walk w as a 2i/- vector of ±l's. It is clear that such a walk has v + I "downturns" 
(vector entries of -1) and v — 1 "upturns" (vector entries of +1). Therefore, specifying the locations of the 
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downturns determines w. From this we see that the set of all such walks, {w} is in one-to-one correspondence 
with all choices of + 1 numbers (the downturn locations), {ji < j2 < • • • < ju+i} from {1, . . . , 2i/}: 

v+1 



Thus (|07|) becomes 
(4.28) 



bl+ei(w)bk+i2{w) ■ ■ 

{w} 



.7,2 



2m + 1. 



where 

(4.29) eM=jr> 

Some explanation is perhaps necessary here. The sum appearing in H4.28I) is taken over all choices of -|- 1 
numbers, and each such choice is denoted by a member w of the set {w}. As described above, each w may 
be interpreted as a walk on the integer lattice, from k+1 to k-1 of length 21/. Since the structure of the 
walks and all associated counting is independent of the value of k, we will from now on, for simplicity, take 
k=0. The quantity im{w) is defined to be the location on the integer lattice after the mth downturn. It is 
now straightforward to compute that, given a walk with downturns at ji, . . . after m downturns, there 

have been jm — rn upward steps, and m downward steps. Since the starting position is -1-1, 14.29|l follows. 
Inserting these reductions into the Toda equations H4.26(l we have 

'u+l v+1 

(4.30) ' ~ ' 
By Proposition ^2 

(4.31) hUO 



{w} 



/c+£„i(u!) + l 



m— 1 



m— 1 



k+lrn{w) 



Sk 



and the differential equation H4.30|l becomes 



= k 



^^log(Z,(ii, 



'Sk)) 



(4.32) 



{w} 



k+e,„(w)+ii^k+e,r,+ 



-i) - n bk+e,„(w)(sk+e,J 



^From (|4.31|) . and the asymptotic expansion satisfied by the partition function appearing therein, one 
sees that for any fixed integer the quantity {k + £)^^bf,^^{sk+t) possesses an asymptotic expansion in even 
inverse powers of (fc + £), which we express in the following form: 



(4.33) 



bl_^_^{sk+e) 



kj:. 

g>o 



Sk 1 




-2s 



Fundamental Asymptotic Principle: Recall that 

i: the asymptotic expansion satisfied by the partition function is differentiable term by term 
ii: the coefficients Zg{u) are analytic functions of •. 

Therefore, the terms appearing in the summation in H4.33|l may be Taylor expanded (note that £/k is 
asymptotically small) and the resulting multiple series may be re-summed. One may convince oneself of this 
as follows. One truncates the expansion (|4.33|) at order k^~'^^ , and observes that the error is 0{k~^~'^^). 
Then one computes Taylor expansions to order A:"^'' of all quantities appearing in the (now finite) sum, and 
observes that the error is 0{k'^^^^). 

Following this principle, wc now introduce the scheme f{s,w): 
(4.34) f{s,w) = Y.zg{skw''~^)w^'^^k-^s. 
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Using this scheme, the differential equation (|4.32|) becomes: 




where we note that (i^ + 1) ^ v'^X ) ^ ^'^ ^ -j^ ^ = c,y as claimed in the statement of Theorem 12 .21 
From the expansion (|4.34|) one may now confirm that the higher order terms in the expansion of the RHS 

have the form described in Theorem 12. 21 The justification for the form of the coefficient c?^'^'' stated there 

requires some further explanation which we now provide. 

Given the explicit determination H4.29|l of the im. in terms of the walks, it is natural to try to calculate c?y 

by expressing this combinatorial quantity in terms of expectations with respect to the induced probability 

distribution on the family of random walks \w) conditioned to begin at +1 and end at -1. To this end we 

observe that given an ordered sequence of numbers m\ < mi < ■ ■ ■ < selected from the set 

and another ordered sequence ii < ■ ■ ■ < ir selected from {1, . . . , 2i/}, 



(4.36) tt {w\ mf downturn of w occurs at position ij, j = 1, . . . , r} 

«1 — 1 \ / ^2 — il ~ 1 \ / ir — ir-l — 1 \ f 2v — ir 

nil ^1 /V rni — mi — 1 / \ mr — m,._i — 1 ) \ v + \ — mr 



These joint probability distributions are fundamental for the count we will now describe. 
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One sees directly from expanding the products in (|4.35|l . that 



distinct 

^ /p(V) P(V) 



mi,.. .,mp(y) 
distinct 



p(v) p(y) 



e {!,..., 1^ + 1} 



EE E 



mi - 1 J \ 1712 - mi - 1 J \ mp(y) - mp{v)-i -1 / + mp(y) 



rni - 1 y \^ m2 - mi - 1 y \^ "7,p(y) - mp(y)_i -1 J \ v + 1- mp^y) 



W {ij - 2mj + 1)1 



where in the last equality we have applied 14.29|l and H4.36|l . 

The formula for dy'^"* given in Theorem 12 . 21 now follows from the next proposition. 



Proposition 4.2. The inner-most summation may be re-expressed in terms of the following combinatorial 
coefficients: 



coeff of x'^ ^'^^in 

(^ -4- T, 

' dx 



\{ U„-2f](l+X,^) + c) ) ■{l+Xiy'-^---{l + Xrr-'^-'~\l+Xr + lf''-'^\,^=, 



E 

mi < • • • < mp(y) 

r 

Y\ {in - 2m„ + c)"" 



n=l 

21 



«1 — 1 \ / ^2 — — 1 A / ir — ir-1 — 1 \ ( 2v — ir 

mi — 1 / V m2 — mi — 1 / \ m^ — m^-i — 1 ) \ u + \ — mr 



Si S2 ^ rrST^r+l I 



To prove this we start with the LHS which, after expanding binomially, becomes 

nU„-2^(i+.,A)+, -E E ••■ E E 

ri=l \ s=l * / / si=0 S2=0 s^=0 Sr+i=0 

T ') ^ ■■("■"::'"') ( ) ^ ■ ■ i— 

ii— 1 12— ii— 1 V — ir-i— 1 2iy — ir r 

= coeff of x''-''+Mn E E ' ' ' E E H " + c - 2 (si + • • • + • 

Si— S2— Sr— s^-i-1— On— 1 

Zl — 1 \ / ^2 — ^1 — 1 A ( V — V-1 ~ 1 \ 2u — ir 

Si J \ S2 J \ Sr J \ Sr^i 

coeff of x''~''+Hn E E ' ' ' E E Yl (in - 2n + c - 2 {si + ■ ■ ■ + s„))' 

Si— S2— Sr—0 s^^i— On— 1 

Sl / V '^2 /V J \ Sr+1 

r 

E Y[{in-2n + c-2{si + --- + Sn)T" ■ 

si-f hSr+i=v-r+l n=l 

ii — 1 A / i2 — H — ^ \ ( ir ~ ir-i ^ ^ \ ( 2v — ir 

Sl J \ S2 J \ Sr J \ Sr+1 

where in the first equahty we use the fact that the monomials in xi, . . . ^Xr+i are each simultaneous eigen- 
functions for all the differential operators — 2 'Y^=i{^ + af~) ~^ ^ ^^'^ therefore the action of the product 
of these operators on each eigenfunction may be replaced by the product of the corresponding eigenvalues 
times this function. 

We next make the following change of variables in the summations 

51 = mi — 1 

52 = 7712 - TOi - 1 



Sfc = mfe - 777fc_i - 1 
Sr = fnr ~ TUr-l — 1 

which then requires that 

Sr+1 = i^ — 7' + l — (si + -- - + Sr) = v — r + I — [rrir ~ r) = ly + 1 — rrir 
and this completes the proof of the proposition. 

4.2. Differential Equations for eg{t) . We will now show that there are differential equations whose 
solutions determine eg{—s) written in terms of the functions Zg{s). We will show that these differential 
equations can be used to derive expressions for the Taylor Coefficients of eg(— s), Hg'^'^ (j), and verify Theorem 

ESI 

Recall the Hirota relation (|4.4(l . Restricted to the level 2i/-time flow this becomes 

rfc+i(e)Vfe-i(0^ 

and is related to the partition function through 



(4.37) 6.(0^ = ^^^^^^^-j^^^b.iOf 



2 Zk{-2^k'^-') 



ZkiO) 
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The logarithm of (|4.37|l produces a hierarchy of difference equations (discrete Hirota relations): Setting 
Sfe = 2k'^~^^ in order to write 



(4.38) 



Zk{0) ' 

and applying a logarithm to both sides of equation (|4.37() one obtains 



(4.39) log 



Zk+i{—Sk+i) 
Zk+iiO) 



2 log 



Zk{0) 



log 



/ Zfc_i(-Sfc-i) 

V Zk-i{0) 



iog(fcfc(e)')-iog(&4o)2). 



We will next describe a continuum limit of these discrete equations which will generate a hierarchy of 
differential equations for the functions eg{t2v) which appear in the asymptotic expansion (|1.7|) . 

If we shift fc 1-^ A: ± 1, then st shifts: Sk ^ Sk±i = Sfe(l ± . With this substitution the left hand 

side of H4.39|l becomes 



(4.40) log('^^±i(^#±i^V21og^^^(-^ 



log 



Zfe_i(-s(l-l/fc)-i) 



Zk+m J ^\Zk{0)) Zk-i{Q) 

Note that from Theorem II. 71 each term in 1)4. 40() possesses an asymptotic expansion of the form 



(4.41) 



log 



Zk+r{-s{l + r/kY-^) 

Zk+r{0) 



9>0 



On)' 



with r = — 1, 0, 1. 

The expression (|4.40|l is a centered difference. We utilize a lemma describing centered difference expan- 
sions: 



Lemma 4.3. If Q is an analytic function, then 



(4.42) 



Gil + z) - 2g{i) + g{i-z)^Yl ^^^'"^(1) 



n>l 



2n) 



as z —^ 0. 

To apply this lemma we note that (|4.4U|) has an asymptotic expansion of the form 



log 



Zk+ii-sii + i/kr-^) 



Zk+iiO) 



2 log 



Zk{--s) 
Zk{0) 



g>0 



20^(1) + gg(i 



log 

1 



Zfc_i(-s(l-l/fc)'--i) 
^fe-i(O) 



where 
(4.43) 

Using Lemma 14.31 we have 



Ggiw) = W^-^^Cgi-w'^-h). 



(4.44) log 



Zk+ii-Sk{l + j:r-')\ ^, fZk{-Sk)\^^ fZk-ii^Skil-^r-'] 



Zk+m 



2 log 



log 



J2 fc'"'""'" 

3>0 n>l 



Zk{o) J ' "° V ^fc-i(o) 



In terms of the coefficients Zg defined in Proposition ^3 the right hand side of H4.39I) becomes 



(4.45) 



log[j2k''<^^gis) 

^a>o 
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We get ode's which determine eg{—s) by equating Hke orders in equation (|4.39|) after making the sub- 
stitutions (|4.44() and H4.45|l . The equations so derived are 

^ (2n + 2)!9u;(^»+2) \w^~'^^-' eg,,,i-nj-^-^ s)] = the fc-^ term in log ( ^ fc-"z„(.) ) . 



n=0 



which establishes the hierarchy 1)2. 9|l . These equations are recursive in the sense that we need to know z„ for 
n < g and Cm iov m < g in order to write a closed differential equation for Cg. We will see how this works 
out in the next section. 

5. Continuum Toda Solutions and Fine Structure of the Asymptotic Partition Function 

In this section we will prove Theorems l2.3l and l2.5l The equations of the continuum Toda Lattice hierarchy 
are given by 1)2.7(1 and (|2.8|l and bk{t) possesses an asymptotic expansion whose leading order, zo{s), is the 
branch of the solution of the polynomial equation 

(5.1) 1 = zo{s) ~ c^szq{s)'' 
which is regular at s = 0. Recall that 

Implicitly differentiating 1(5.1(1 we have the basic identity 

(5.2) ds = t^S-_^dzo. 
which will be used often in following sections. 

5.1. Integrating the Continuum Toda Lattice Hierarchy. Let fg{s,w) — w^^'^^ Zg{w^~^s). The k~^^- 
order equation in 1(2.7(1 is linear in Zg{s): 

(5.3) z'g{s) = (fofgw + i^fo^^fgfotu) \iu=i + Forcingg|u,=i 

where Forcing^ is given by ((2.8() and the entire RHS depends only on Zj,j < g and their derivatives. (Note: 
in what follows, if the context makes it clear that the forcing terms are being evaluated at w = 1 we will 
simply write Forcing^ rather than Forcing^j^^i.) 
Expanding out the derivatives in equation ((5.3(1 . 



(5.4) z'g = Cn ((1 - 2g)z^Zg + (j/ - l)sz^z'g + vz^Zg + v{y - l)sz^ ^Zgz'^ + Forcing^. 
Solving for z'g. 

(5.5) (1 - c,(v - \)sz'(^z'g = ((j/ + 1 - 2g)zl + v(v - \)szl-^z'^) Zg + Forcing^. 

We may inductively assume that Zji^s) is a function of zo(s) for j < g, and as a result Forcing^ is a function 
of zq. One can convert equation ((5.5(1 to a differential equation for Zg as a function of zq using equations 
(Umi and (lOl : 

(5 g) d£g ^ (z.(2-2g) + (t/-l)(2ff-l)zo) _ , Forcing^ 



dzQ zo{iy - {ly - l)zo) ^ c^Zq^^ 
Equation ((5.6(1 is linear in Zg. Let 

Jy U{iy - [l' - 1)U) 

(5.7) = / ((i/-l)-i(i/-(i/-l)u)-i + 2(l-5)u-i)du 



U=Zo 

(5.8) =-log(i.-(i.-l)^/) + 2(l-5)log(u) 

u=y 
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Integrating equation (|5.4|l one has 

Zg{s)= / cxp(G'(zo,?;)) —^Forcmgg{y)dy 

Ji c^y 

2(l-g) i-zo _(-,,_ 1 ^,, 

(5-9) = TS- / r+3- Forcing^(y)dy, 

where zq = zo{s). This estabhshes the basic formula stated in Theorem |^2t3). 

5.2. Singularities of Zg. In this section we complete the proof of Theorem 12.31 Expression (|5.9|l can be 

exploited to determine the singularity structure of Zg{s). We differentiate equation 1)5.1(1 and solve for z'q{s): 

(5.10) = , 7°^' 

(i^- [V- 1)Zq) 

With equation l|5.1U|l we can re-express successive derivatives of zq as functions of just zq by using the chain 
rule. 

Doing so, we find the following result: 

Lemma 5.1. The derivatives of zq{s) have the following form 

■j ^''J+i 
- - l)zo)2j"- 

Equation L5.11]) together with 1^5. ^) give the following result 

(5.12) s^z\^\s) = ° ° — yT—T ■ ( Polynomial of degree j - 1) . 

(v- [v- l)zo) 

The terms in Forcing^ 1(2. 8|l are each of the form 

(5-13) (/fel)u,("l)(/fe2)u.<"2) ■ • ■ (./fc„+l)^('V + l) 



(5.11) 4''^('5) = . _(^''^\\ — ylf^ ' ^ P(^^V™'^^°-^ degree j - 1) 



w—l 

where = w^~'^^ zi^{w'^~^ s) , and rti +n2 + • • • + = 2{g~ki^k2 — ■ ■ ■ — fcjy+i). Expanding the derivatives 
in H5.13|l one finds that Forcing^ is in fact comprised of terms of the form: 

(5.14) s^z^^hl^:\..zl^X^\ 
where 

j = ji + 32^ 

< j < 2(5-fci -fc2 K+i). 

Lemma 5.2. The function Zg^\s) (j derivatives of Zg{s)) is of class iir with singularities occurring only at 

Zq = or Zq — V j (v ~ X) . 

We will first prove by induction that Zg{s) satisfies Lemma 1^?^ see section 1^31 for the initial step. Assume 
that z^"*' (s) for k < g satisfies Lemma 15.21 this assumption, Lemma 15.11 and formula (|5.1|) show that the 
integrand in formula (|5.9|l is a polynomial in functions of class iir with singularities restricted to and 
j^/(i^ — 1) and having rational coefficients with at worst poles located at zq = 0. Thus the integral remains 
of class iir with the stated singularities and H5.9|) is then this integral multiplied by a rational function with 
pole only at — 1). 

We next prove, by induction on j, that the higher derivatives, Zg^\s), satisfy Lemma 15.21 which will 
complete the original induction. A derivative with respect to zq will raise the order of a singularity but will 
not introduce new singularities. Compute 

by the induction step (in j) and formula H5.11|l this expression is of class iir with singularities only at zq = 
and Zq — ^/(^ ~ !)• 
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With the forthcoming verification of the initial step of the induction (subsection , this completes the 
proof of Theorem l2.3l 



5.3. Example: g — I . Recall that fo{s,w) ~ wzo{'w'' ^s) and fi{s,w) — w ^zi{'w'' ^s). The forcing for 
zi is given exclusively in terms of zq as 

(5.15) c;iForcingi = c-'F^O] = ^/o^/o™ + ^^^^/o^" Vo./o.. + ^^^^^^^|^/o "Vo^ 
The entire expression is evaluated at w = 1. We find that 

(5.16) c-^FmU. - K^m^^^^r ^ ^(9.^ _ 17. + 6).z„^4 + '-^^^s^44 

Inserting (|5.16() into equation (|5.9|l with g — 1, (writing all derivatives out as rational expressions of zq) 
yields 

,^ „N \)v[zq - l)zo(-J^2 - 2zo + VZq + v'^Zq) 
lo-j-'j ^il'j — 7777 7 T^; — u ■ 

^From this we see that not only does z\ it) satisfy Lemma 15.21 but in fact this expression is a rational 
function, without any of the potential complexities of the iir class. In fact all the examples we have worked 
out here for low values of g yield rational expressions in zq for Zg. Moreover, the poles of these expressions 
are restricted to zq = vjiy — 1). 



Remark The rationality of H5.17|l seems to depend very sensitively on the precise coefficients appearing 
in Forcing^. We illustrate this by looking at the contributions that particular terms of H5.16|l give to H5.17|l . 
The integrand in integral H5.9|) has a term of the form 



(z/- (j/- 1)Z0) 2 .-1. /^2 

-S Zq 1^07 



(^0 - 1)^ 



[y - [y - l)zo) 



where we have used equation (|5.10|) to eliminate Zq. This term produces a logarithmic singularity at v jiy— 1); 
therefore one of the other terms must cancel this one to eliminate the logarithm. □ 

^From the representation H5.17|l we can compute the Taylor Coefficients of zxify using the technique 
outlined in formulas H3.44I3.46|I : 



.(1) ^ ^ 
27ri 



[y-[y- l)u)u^J-i 



zi(zo u)du 



du 



(5.18) 



2^1 l2{u-iy{v-{v-l)uY 

2TTi 



12(1 -{ly- l)w)^wi 



dw. 



In the last line we have made the change of variables w = u—1. With the integral representation H5.18(l one 
may conclude that 



Taylor Coefhcient of 



{{v ~ ^2) + {v + 2)w) {w + ly^ 



(1 



l)wf 
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5.4. Example: g — 2. The forcing terms for Z2 (|2.8|) are: 

(1) The first term is the summation in (|2.8() : 

(2) The middle term is f[''\2]: 

^1 ''[2] — "^/o flfowivw + ■^/o/luiioto ^ ^ /o flfowfoww ^ ^ /o flwfoww 

+ 3/0 JOwJlww ' 22 /l/Oui ^ ^ •'0 /Ouj/liD- 

(3) The last term arises from F2''''[0]: 

"-^2 i^JL=i— J^^^ Jo/Otu(5)H ^^^^ Jo /0t„/0u,(4) H Jo JOwwJOwww 

t/(z.-l)(:.-2)(23^-20) , , t/(^-l)(^-2)(62^-65) , . 

360 J ^wJ Owww ' 720 J Ow J Oww 

K'^-l)(;.-2)(z.-3)(16z.-19) 3 3 
360 J J OwwJ 

v{v - - 2){v - 3)(i/ - A){hv - 7) ^,„4 ^ 5 
+ 1440 

The three terms make up the forcing terms for Z2; each of them is expressed as a function of z^. Substi- 
tuting this explicit representation into equation (|5.9|) with g = 2 we find that 

^2{t) ^-^{^ - - 1)^0 [(21^^ - Ui^' + 24i.«) 

+ {-l2v^ + 148t/'' - 546i/^ + 758i/^ - 2b2v^ ~ %Qv^)zo 

+ (264iy2 _ _^ 255511^'' - 500i/^ - 17891/*^ + 840t/^ + UAv^)zl 

+ (-536i/ + 13961/2 + 912i/3 - 4596;/^ + 2A92v^ + 12961/*^ - 868i/^ - %Qv^)zl 

+ (168 + 234i/ - 1467t/2 + 558:/^ + 1902:/^ - 1446t/^ - 267:/*^ + 294i/^ + 24i/^)z;5] 

•(i/-(i/-l)zo)^'. 



Let 



(i/-2)(i/-3)(j/-4)(5j/-7) 

+ 2(iy - l)(i/ - 2)(168 + Mv - 246J/2 + 7?,v^)w 
+ {v- lf{v - 2)(-504 - llb^v + 288i/2 + 'i<dlv^)w'^ 
+ A{v - 1)^(-168 - 604i/ - 190i/2 + 2%^v^ + 77j/'')u;3 
+ 3(i/ - 1)^^(56 + ZQ2v + 383Z/2 _,_ i30j,3 gj^4^^4 



We compute the Taylor Coefficients of 22 (i) with a loop integral: 

>(2) _ Ci / 1)U)U-^-^ 



=tlL, (U-1). + 1 .2(.0=.)rf. 

= The w^~^ coefficient of r2(w). 

5.5. Example: g — 3. To compute Z3 we compute the forcing as above. Here we outline how the terms 
contributing to this forcing are found and then present the resulting explicit expression for Z3 . There are six 
terms in formula H2.8|l when .g = 3: 
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(1-2) The first two terms occur in tlie summation in H2.8|l when g = 3. One term contains three /I's and 
is 



g Jo /l /Oto + 2 •'0 /l/liu- 

The other, containing an /i and an /2, is 

(3-4) The next two terms are contributed by F^'-'^-' [4]: this expression represents the fourth order contribu- 
tions of f['^\ The terms involving two /I's are 



/o /l fOwww + 



12 



/o fiwflww + 



/o flfOwf- 



W J J- WW 



WW ' 



^.(i.- 1)2(^.-2) ^,_3 



/o /l /Oto/o 



K;^-l)(;.-2) 



/o /Otu/l 



K^^-l)(^.-2)2 „ 



W ' 



K^-l)(i/- 2)2(^-3) 4 2 3 

2^ J a JiJOw 



The terms containing one /2 are 



'Jo Jf^wJ2ww o •'0 J2wJ0t, 



^fof 2www ~^ 5 ./'2y Ou''(i)tjj ~^ 2 



^{^~^)%u-2f f f ^(^-l)(^-2) ^,_,^2 , , ^(^-l)(^-2)% ,_3^ ^3 
Jo J2j0wJ0ww ^ ^ Jo J0wJ2w ^ 77; /o /2/0^. 



3 



12 



(5) The fifth term is 



K2^-l) i.2(2z.-l) , K^-l)(18^-ll) 1 
-^^2 l^J — 7777; /o/ltu(-'>)H 7777 Jo 7l/0u,(5) H 5^7; Jo JOwJlwW 



120 



120 



360 



+ 



+ 



'-l)(18z.-ll)^,_i 



360 



/o fiwh 



K^^- 1)2(18^.-11)^^ 



z/-l)(6z.-5) 1 

J OwwJ ^ wwiv 



360 

j.(^_l)(6^._5) 



/o flfowfo 



72 

i.- 1)2(61. -5) 2 

— Jo JiJo 



72 



/o /l luui /o 



72 



^(^-l)(^-2)(23^-20) 2.2 r 

Jo J^wJ^www 



j.-l)(^.-2)(23^._20) 2, , , 

Jo JOwJlwJO 



360 

i.(;.-l)(i.- 2)2(23^.-20) ^2 



180 

^-l)(^-2)(62z.-65) 2. . . 

^0 J OtoJ OwwJ ^ww 



360 



"/o /i/otu/o 



360 

z.-l)(z.- 2)2(621. -65) ^,_3 



^(^-l)(^-2)(62^-65) 2. .2 

^^20 '''^ J ^w J Oww 



720 



,2 , ^(^ - 1)(^ - 2)(i. - 3)(16z. - 19) ^,_3 ^ 3 ^ 
Jo JiJOwJOww ^ 5^7; Jo Jo wJ ^WW 



360 



^-l)(^-2)(^-3)(16t.-19) 3.2 . f 

-j^2Q •'0 J OwJ ^wJ Oww 

z.-l)(i.-2)(z.- 3)2(16^.-19) ^,_4^ , 



360 



/l^ — 4 prop 
JlJOwJO 



WJ ^ww 



^-l)(^-2)(t.-4 )(5z.-7) ,._4,4 , , t.(^-l)(^-2)(^-4)2(5;.-7) ^,_5^ ^5 
Jo JOwJiw "1 TTiTT; Jo JiJOu, 



288 



1440 
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(6) The last term is 
^(''^[O] = — C ' ^--^^ 



7!2 






2) 


7!6 






2) 


7!3 




v{v — — 


2) 



(105 - 210;.+ 112z.2)/--i/o^^/o^<,, + ^^(^1^(539 - 1050;. + 504z.2)/--i/„^^,3,/o„(4, 



7!2 

^-l)(^-2) 



7!6 

^- l)(j/-2)(i/-3) 



7!6 

V - l){u - 2){v - i) 



7!2 

i/-l)(i^-2)(i/-3) 



7!6 

(86 - 163i/ + 86z.2)/--2/o^/o^(5) 

(1387 - 2416J. + 1044z.2)/--2/o^/„^^/q^(,, 

(467 - 803^. + 327z.2)/--2/o^/o^(3, 
(1456 - 2359^. + 882i/2)/--2/oL/o^(3) 

(410 - + 270z.2)/--3/o' /o^(4) 
(857 - 1309J. + 458^.2)^.-3^^2^^^^^^^^^^^ 
(1327 - 1932J. + 620i.2)^.-3^^^^^3^^ 

11211. + 3591.2)^.-4^^4^^^^^^ 



7!6 

!.-l)(i.-2)(i.-3)(i.-4) 



7!12 

V - - 2){v - 2,){v - A) 



7!12 



(2431 - 33151. + 974i.2)/o--VoL/oL 



V - l){v - 2)(i. -'i){v- A){v - 5) 



(228- 290i. + 77i.2)/o^-Vol/o 



7!6 



V - l){v - 2){v - 3)(i. - A){v -b){v- 6) 



(124 - 147t. + 35i.2)/o^-Vo' 



7!72 



The entire forcing expression is a rational fmiction of zq with singularities only at zq — vjiv — 1). We 
insert these forcing terms into equation H5.9|) with g = 3 and find that: 



^3(0 



_ v{y-\) zo(zo-l) 

362880 {v- (i.- l)zo)i4 
[(;. - 2){v - 3)(i. - 4)(i. - 5)(i. - 6)(124 - 147t. + 351.^) 
+ (i. - 3)(i. - 2)(z. - 1)(104160 + 475841. - 3325501.^ + 270697;.^ - 83226J.'* + 8923j.^)(zo - 1) 
+ Z{v -2){y- 1)2(312480 + 744980;. - 12457501.^ + 373091^^ + 1085920^^^ - 485414z.^ 

+ 67225z.^)(zo - 1)^ 

+ (i. - 2)(y - 1)3(-1562400 - 7251840;. + 290690^^ ^ 11468057:.^ - 28240781.^^ - 3154302;.^ 

+ 10786631^6) (00 - 1)^ 

+ (z. - 2)(i. - 1)'*(1562400 + 10781280^ + 12588010;.^ - 10677353i^^ - 11255921^'' + 3006363;.^ 

+ 1779986^6) (zo - l)'^ 

+ 3(i. - 1)^^(624960 + 5411808Z. + 10100796z.2 - 1315908^^ - 93716951^^ - 973573;.^ + 1835799:.^ 
+ 308858^^) (zo - 1)^ 

+ (j. - l)6(-624960 - 6823584i. - 20098900^^ _ i6851720i.^ + 3867117!.^^ + 83564421.^ 

+ 2223760J/6 + 119824z.^)(zo - 1)*^ 
+ h{v - 1)^(17856 + 235296J. + 939236;.^ + 15050641.^ + 1032603;.'* + 285860^^ + 24472;.^ 
+64^^)(2o-l)'] . 
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Let 

cji^jv- 1) 
re, u; = (w + I) ■' 



{ly - 2){iy - 3){iy - 4)(j/ - 5)(z/ - 6)(124 - 147j/ + SS;/^) 



362880 

+ (i^ - 3) (i^ - 2) (i/ - 1) (104160 + 475841^ - 332550t/2 + 270697i^3 _ 83226;^'^ + 89231/^) w 
+ - 2)(iy - 1)2(312480 + 744980i^ - 1245750i/2 + 373091;^^ + 1085920;^'* - 485414i/''^ + 67225i^^)w2 
+ {v-2){iy~ 1)3(-1562400 - 7251840;/+ 290690i/2 + 11468057z/^ - 2824078i^'' - 3154302i/^ 

+ 1078663z/<')w3 

+ {iy-2){u- 1)''(1562400 + 10781280i/ + 12588010i/2 - 10677353^^ - 11255921i^'' + 3006363;/^ 

+ 1779986z/6)«7'' 

+ 3{v - 1)5(624960 + 5411808i/ + 10100796^^ - 1315908;/^ - 9371695^'' - 973573^^ + 1835799;/^ 
+ 308858i/^)m;'^ 

+ (^ - 1)^ (-624960 - 6823584Z/ - 20098900:/^ - 16851720;/^ + 3867117:/^ + 8356442:/^ 

+ 2223760^^ + 119824^^)u;6 
+ 5(^ - 1)^(17856 + 235296;/ + 939236;/^ + 15050641/^ + 1032603;/'* + 285860;/^ + 24472;/^ 



+ 64i^^)u;'^ 



(l-(;/-l)u;) 



-13 



We compute the Taylor Coefficients of 23 (t) with a loop integral: 



A3) ^ 

2m 



The w^~^ coefficient of r^{'w) 



•23(^0 = u)du 



5.6. Determining eg{—S2u)- Theorem 12 .41 gives a second order differential equation determining eg{~S2u)'- 



(5.19) 



)] = dr 



ivers„ 



where the entire equation is evaluated at w = 1 and where the drivers are the right hand side of equation 



E 



a2n+2 



■'^ (2n + 2)! dw^''+^ 



w 



2-2{g-n) 



Cg-ni-w" s) 



+ the k ^3 term of log 



J2 T2^^n{s) 



Expanding the two derivatives on the left hand side of (|5.19f) and setting w = 1 we arrive at the basic ode 



(5.20) (2 - 2g)(l - 2g)eg{-s) - (j. - l)(i. + 2 - Ag)se'g{~s) + {jy ~ l)^s^e'^{-s) = drivers^. 

If g 7^ 1 one solves this equation by integrating factors as follows: Multiply both sides of (|5.20|) by s^'^ 
with 



71 = 7i = 



(+) _ 3 - ;/ - 2ff 



i/ - 1 



or 7i = 



(-) _ 2-iy-2g 



V - 1 



(5.21) 
to arrive at 

(5.22) (2 - 2g){\ - 2g)s'''eg{-s) - {v - + 2 - Ag)s'<^+'^e'g{-s) + [v - if s''^+'^ e'g{- s) = s'^MriverSg 
Let 

A= (2-2g)(l-2g)/(7i + l). 
We then integrate equation (|5.22|l once, 



(5.23) 



As^^+^eg{~s) ~{y- \f s''^^'^ „{- s) = / si'^idriverSg(si)dsi + K{ 
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where K[ is a constant of integration. Multiplying by s"^^, 

(5.24) As'''+''^+^eg{s) - {v - if s''^+''^+'^e'g{~s) = s'^'^ / si''MriverSg(si)dsi + s''^K[, 

Jo 

we see that if 72 satisfies 

(.^r. ^ (2-2g)(l-2g)-(^-l)^(7i + l)(7i+2) 

(71 + ly 

one may solve equation H5.24|) for eg{s): 

(5.26) eg{-s) = L__s-71-72-2 / ^72 / e7lrl,.hrPrB_reoW«oW«, 4- ft', Q-71-1 _|_ ft'^e-71-72-2 



/ / S2^driveTSg{s2)ds2dsi + Kis ^ + K2S 
Jo Jo 



where Ki and K2 are constants of integration. 
Switching the order of integration in (|5.26(l 

(5 27) e„(-s) = ^ g-7l-72-2 f f a72 ^7lHrW70rc /'o„Wo,^o„ J- ft-, 0-71-1 _i_ R'„a-7l-72-2 

one can compute the si integral in H5.27|l 

1 1 



/ / sj^ S2^dnverSg{s2)dsids2 + Kis '^^ ^ + K2S 

Jo Js2 



(5.28) e,(-,s) 



- 1)2 72 + 1 
^ Kis-'''-^ + K2S 



'f'^^ f s]MriverSgds2-s"'''"^'"^ / s]^+^=+^driverSgds2 
Jo Jo 



s 

-71 — ("2-2 



Applying an integration by parts to the remaining integrals (integrating the power of S2 and differentiating 
the drivers (with respect to S2)): 



es(-s) 



1 1 



(^-1)^ (72 + 1) 



-71 — r2— 2 ps 



(5.29) ' ' I sI'+'''+'(driverSg)'z^ds. 



71 + 1 71+72 + 2 



drivers^ H / S2^^ (driverSg)*ZgC?S2 

71 + 1 Jo 

+ Kis-'^^-^ + K2S-'^^-'^^-'^. 



71 + 72 + 2 7o 

We may now change variables and integrate with respect to y = ^0(52) in H5.29|) 

72 + 1 



egi-s) 



(;.- 1)2(72 + 1) 



driverSo 

(7i + l)(7i+72 + 2) " 



1 f r \ 71+72+2 „zo / -|X 71+72+2 



71 +72 + 2 V^o - 1/ 7i V Ci/J/'' 

71-72-2 



(driverSg)*(i?/ 



(5.31) +Xis-^i-i +i^2S 

where K\ and -fC2 are constants determined in either of the following ways: (1) by the requirement that Cg 
be a locally analytic function of s or (2) by the evaluation of Cg for low values of v through its combinatorial 
characterization. 

Direct calculation shows that either choice of 71 in (|5.21|) produces formula (|2.1U|) in Theorem 12. 51 
5.7. Example g = 0. When g = we have 

(-) ^^-2 

71 =71 ' ^ 

V — \ 

then 

v-1 



The driver when 5 = is just log(2o)- The two integrals in l|5.31|l may be evaluated separately: 



(y-l)V(- 



l)y-(2l'-l)/(l'- 



'^dy= / (1 



-1x1/(1.-1), -2 



y dy 



(5.32) 
and 



_ i)2/(.-i)j^-(3.-i)/(.-i)rfy = / (1 - y-'f/^-^-'-'y-'dy 



(5.33) 



:|^(l-Z-)(''+l)/(^-).- 

('^ + 1) Jl 

■(^ + 1) 



2u/{i^-l) 



Plugging ifO^ and into l(OT|) with g = yields: 

eo(-.) 4 log(.o) - ^(.0 - 1) + i^(.o - 1) + l^^^i^o - 1)^ 
I1 / X ,x ^ 3(1. + 1) 

where Ki and are both by the analyticity condition. This formula for eg agrees with the expression 
for eo(s) found using the equilibrium measure. 



5.8. Example: g — 1 . The drivers for ei are 

Zi 1 



zo 12 9w4 

evaluated at u; = 1. As with the Forcing terms for the Zg this formula can be expressed as a rational function 
of Zo with poles at zq = and zq = v jiy — 1). 

We find that 5 = 1 is the exceptional case for the method used in formulas H5.22ll - (|5.31|l . We compute it 
directly from equation (|5.20|l with g = 1. The integrating factors are (71,72) = (—1, —v/{v — 1)), differing 
from the general choice of integrating factors (|5.21|l and H5.25|l . 

We find the integral formula: 



ei(-s) 



1 



(^-1)2 
1 

1 
1 



S2"Mriversi(s2)ds2dsi +Ki+ i^as^^'""^^ 



/o Jo 

PS ps 



S2"Mriversi(s2)dsids2 + Ki + K2S^/'-''-^'> 



JS2 



/ 82'^^'"'^ ^■'driversi(is2 — / S2 ^driversi(is2 
"'0 "'0 



Zois) - 1 f^ois) ^ l)y) . 

' III — driver Si ay 



r driversi ay 

y[y-^) 



c^y"-^' 



= - ^ log (i/ + (1 - i.)zo) + Ki + K^s^'^"-^^; 

where we have switched the order of integration, computed the si integral, and changed coordinates to 
integrals with respect to y = zo(s2)- The final step is a direct computation and is left to the reader. 
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When t = 0, log (ZAr(t)/Zjv(0)) = 0, therefore ei(0) = and so Ki = 0. We notice that if z/ 7^ 2 then 
K2 = by the analyticity condition, otherwise we will have to determine A'2. The computation we will carry 
out in Section [5 . 1 01 shows that when v = 2, 

ei(-s) = s + 0{s^). 
Our expression for ei(— s) has the expansion (when v ^ 2) 

log [2 - zo{s)] + K2S = (1 + K2)s + 0{s^), 
therefore K2 = 0. This proves formula (|2.11ll . 
5.9. Example: g — 2. The drivers for 62(5) are 

evaluated at w = 1. As with the Forcing terms for the Zg and the drivers for ei this expression can be 
expressed as a rational function of zq with poles at 20 = and zq — v / {v — 1). 
For g — 2 wc find that 

(+) + 1 

and 

V 

72 = - 



We compute the integrals in equation 15.31|l . 

e2(-s) =2^(^ - l)(^o -\)(y-[v- 1)20)"' ((-^' + hv' + 8^5) 

+ (-1^2 + 41z^^ - 24z^'' - \<oU^)zq + (44zy - SQt/^ + 54iy3 - iri/"^ + ^v^)zl 

+ (-12 - 12;^ + 108i^^ - 132zy3 + A&v%l + (-12 + 48;^ - 12v^ + 48;^^ - \2v%^ 

We notice that K\ and if2 will be zero by the analyticity condition unless = 2, 3, or 4 in which case they 
will have to be determined by some other means such as that illustrated in the next subsection. 

5.10. Evaluating the constants of integration . We will now outline a method for determining the 
constants of integration in formulas H2.10|l and (|2.11|) of Theorem 12. 51 

We know, from |15.j , that eg{—s) is analytic in s in a neighborhood of s = 0, therefore we notice that 
for many values of 71 both K\ and Ki will vanish to preserve analyticity. However there are some values 
for which they will not: Values of v and g which conspire to make [2g — 2)/(z^ — 1) or {2g — l)/{v — 1) be 
positive integers. To find these constants, when necessary, one may rely on the combinatorial interpretation 
of the Taylor coefficients of Cg described in Theorem ll.3l Here we outline the method that we have used to 
evaluate these constants for 5 < 3. 

The set of maps D having vertex set Kq{D) of fixed cardinality at specified vertices can be placed in 1-1 
correspondence with a class of subgroups of the permutation group Sd where d, twice the number of edges, 
is determined by the cardinality and valences of Kq{D) 0]. For simplicity and also because it is the case of 
relevance for us, we take Kq{D) to consist of n vertices each of valence 2v. For this class d = 2vn which 
is the cardinality of the set of darts of K{D). A dart is an oriented edge. Equivalently, we may define the 
abstract set of darts associated to this class of maps as 

= U fi2, 

where 

rJi = {(w, e) : w is a vertex and e is an edge with two distinct vertices one of them being v} , 

and 

^2 = {{v, e, ±) : u is a vertex and e is an edge with a single vertex u} . 
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The element (w,e) € fli represents the dart based at v and going along e. The element (?;,e,±) S ^2 
represents the dart based at v going along e in the counterclockwise (resp. clockwise) orientation. For each 
edge there are two darts, therefore = d and we can think of Sd as acting by permutations on the set fl. 

Given a map, {K{D), [l]) with Kq{D) specified as above, we define a subgroup of Sd generated by two 
permutations {cr,T). The orientation on X induces (via [l]) a cyclic ordering on the darts attatched to each 
vertex; the first permutation a is given by this action. Explicitly, a maps the element {v, e) S Vli to the 
element (u, e) € Jli or (d, e, ±) € ^2 where e is the edge counter clockwise in the orientation at v from 
e. Likewise a maps the element (w,e,±) S O2 to the element (w,e) € Oi or (w,e, ±) G ^2 where e is the 
edge counter clockwise in the orientation at v from e. The second permutation r, is given explicitly as the 
permutation which acts on f2i by sending [v, e) to (ti, e) where v is the other endpoint of e; r acts on VL2 by 
sending {v, e, ±) to {v, e, =f). 

^From these descriptions one sees that: cr is a product of disjoint 2j/-cycles, with each cycle corresponding 
to a unique vertex in Kq(D)] and that t is a product of disjoint 2-cycles, with each 2-cycle corresponding to 
a unique edge in Ki{D). 

Conversely, given a subgroup presented as above which also has the property that the group acts transi- 
tively on fi, one may construct a unique map. The transitivity condition insures that the underlying map is 
connected. The permutation r determines how the vertices are connected through edges to define the graph. 
The permutation a gives the orientation of the edges about each vertex. Together these two permutations 
determine the surface X and the embedding class [l] . 

The punchline is that what we now have is an algorithm for computing the coefficient of s"/n! in eg{—s) 
for finite g, v, and n: Let d = 2vn. Fix cr to be a permutation formed by a disjoint product of n 2^'-cycles in 
Sd- Then we choose each disjoint product of vn 2-cycles, r in 5^. Check if (cr, r) is connected (by verifying 
that the orbit of (cr, r) • 1 is all d letters). If (cr, r) is connected compute the genus by Euler's formula 
(1 — v)n + F = X = 2 — 2g, where F (the number of faces) is given by the number of cycles in cr o r. The 
details of this calculation together with some examples and a generalization to the case of unoriented maps 
are in a forthcoming paper '23^. 

The algorithm gives the following values for the coefficient of s-'/j! in eg{—s): 



9 = 


1 






i 


= 





9 = 


1 


V = 


2 


3 


= 1 


1 


9 = 


2 


V — 


2 


i 


= 2 





9 = 


2 


V — 


2 


3 


= 3 


1440 


9 = 


2 


V — 


3 


3 


= 1 





9 = 


2 


V = 


4 


3 


= 1 


21 


9 = 


3 


V — 


2 


3 


= 4 





9 = 


3 


V — 


2 


3 


= 5 


58060800 


9 = 


3 


V — 


3 


3 


= 2 





9 = 


3 


V — 


5 


3 


= 1 





9 = 


3 


V = 


6 


3 


= 1 


1485 



Then Ki and K2 are chosen so that the j'th coefficient of eg{—s) matches these numbers. 



5.11. Multiple times . The constructions carried out in this paper extend to multiple even time parameters; 
this is a reffection of the commutativity of the underlying flows in the Toda Lattice hierarchy. However, the 
expressions found are more complicated and are less easily reduced to closed form than the monic time cases 
we have considered thus far. 

For brevity we will show how our results extend to two times: t2vi and t2u2 1 where vi and 1^2 are positive 
integers. Theorem 12 .21 becomes 

Theorem 5.3. In the limit as k ^ 00, bk{S,2v 1,^,2^2)^ ^'J* a valid asymptotic expansion of the form 

bk{£,2wi,^2v-2f - k{zQ{s2vt,S2u2) + Zl{s2ux , 82^2) + -^Z2{s2ui , 82^2) + ■ ■ • ) 
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where 821,^ — 2k'^' • The terms of this expansion are determined by the following partial differential 

scheme: 



evaluated at w — 1 



where 



f{s2ui,S2u2]w) = fo{s2u^,S2v2^lW) + fi{s2ui, 82^2'- + (s2iyi , S21/2 ; + ■ • ■ , and fg has the form 

fg{s2„^,S2u2;w) = w'^~'^^ Zg{w''^'^ S2i,^,w''^~^ S2i,^). 

This partial differential scheme yields a hierarchy of partial differential equations in the same manner as in 
Theorem \2.^ The functionals Fg'^^'^ are identical to the ones found in Theorem[ 



The Zg(s2i/i, S21/2) piety the role of auxiUary functions in computing eg(—S2i/i, —521^2)- The differential 
equation determining eg(— S2j/i, — S21/2) remains largely unchanged from H5.19fl because the construction in 
section applies for the multi-time case as well: 

Theorem 5.4. The g'th equation in the hierarchy of equations governing eg{t2ui,t2u2) 
(5.34) 



[w^-^^eg{-w^^'^S2., , -^^-'82.,)] 



Q2n+2 



^ (2n + 2)! 9u;2"+2 

n— 1 ^ ' 

the k~^^ term of log 



„2-2(g-«) 



eg^ni-w"^ 821,1, -w""^ 



00 ^ 



Equation (|5.34|) determines eg{—82i,-^, —821,2) from a second order partial differential equation of Cg with 
forcing terms depending on e„, n < 5, Zn, n < g, and their derivatives with respect to 821,^ and 521^2- 

Theorem 5.5. The RHS of will henceforth be denoted by driverSg(s2i/i , ^2^2) ■ The solution of J^^j) 

may be represented as 

(5.35) eg{-82i,i,-82i,2) = (vi - 



-(1-23)7(1.1-1) 



'21^1 



S2„l 



-(2-2g)/(i.i-l) 



= 2i/i 



(S21.1) 



(2-i'i-2g)/(i.i-l) 



driversg(s2i/i,us2^^ ^^^^"^ ^^)ds2i. 



(3-i/i-23)/(i/i-l) 







J • /" '(^^2 — 1)/(!^1 — 1)\ I- 

driversg(s2,.i,us2i.^ )as2i- 



+i^i(.)4i:'^^/<^^-^^+i^2(«)4t; 



(l-23)/(z.i-l) 



where u = 821^2^2!^'^^ ^^^f where Ki{u) and K2{u) are analytic functions of u in a neighborhood of 

w = 0. 

To prove Theorem 15.51 first we will show that 2:0(521/1,521^2) is given implicitly as the solution of an 
algebraic equation with coefficients depending on S21/1 and S2u2- Then we demonstrate that Zg{82vn 82^2) 
are functions of 2:9(521/1 , S21/2) and zo(s2i/i , 0). Finally we compute eg(— S21/1 , —82^2) by integrating the partial 
differential equation in Theorem l5.34l 

The first order terms in the hierarchy of Theorem 15.31 are the pair of equations 



dS2yi 

dzp 

dS2i/2 



c^iZqI yzQ + {vi - 1)82 
c^jZq^ ( zo + {vi - 1)82 



dzQ 

rfs2i/i 

dzQ 



{V2 - 1)S2V2 , 

as 



dzo 



ds 



2vi 



{V2 - l)S2u2-r 

as 



2V2 

dzo 



2 1/2 



or in vector notation 
(5.36) 



M 



( dzQ 
dS2v-l 
dzQ 
dS2„2 . 



,1/2 + 1 I ' 
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where 



Invert M in (|5.36|) to find a pair of ordinary differential equations for zq: 



dS2v, -1 + {Vl - l)s2,yiZo' + Cy2(^2 - l)s2,.2Zo^ ' 

wfiich, with initial condition za{Q, 0) = 1, is satisfied implicitly by the solution of 

(5.37) 1 = ZO - C^,S2i/i^o' - C1.2 52^2^0^ 

which is regular at (s2i/i, 521^2) = (0: 0)- 

The (7^^ pair of equations in the hierarchy of partial differential equations in Theorem 15.31 is 

/ dz 

(5.38) m(^- ) =z,G, + F„ 




where 



G 



' " U.2(^2 + 1 - 2g)z^,' + C,,V2{yi - l>2.,Zo'''^ + C,,V2{U2 - l)s2uA''^M- 



^ ^ f Forcing^'^i) 
^ l^Forcing^''^'^ 

and where Forcing^'^''' is given by (|2.8|l with 

fg{s2„i,S2u2;w) = W^^^'^ Zg{w''^^'^S2^l,w''^^^S2:,2)■ 
lt is useful at this stage to make the following change of variables: let zo(s2iyi,0) = i/q and denote 
1^1,321/2) zq. The constraint equation H5.37I) becomes: 

(5.39) 1 = yo - Ci.iS2i.i2/o' 

(5.40) 1 = Zo- S2ui z'^'' ~ S21.2 Zq" , 

and one can solve equations H5.39|l and (|5.40() for (521/1,521^2) functions of (?;o,2o)j 

(5.41) S2u^^{vo-l)Mo'r^ 

(5.42) 52.2 = (2/0^(^0 - 1) - ^o^(yo - l))(c,2yo'^;^o'')-\ 

and then differentiate equation (15.391 and 15.401) with respect to S2;y, and 82^2 £^nd solve for f^" , f-^° and 
as functions of (yo, 20): 



dS2ux 

dyo 

ds2u2 

dzp 

dS2vi 

dzt 



1'iCi,iS2i^iZq^ ^ - V2Cu232u2Zq' ^) ^ 



— ~ " " - 1^101.1521.120 - i^2Ciy2 521/2^0 ) ■ 



d52. 



-^1/2 "^0 



Changing variables in the system of differential equations H5.38|l to differential equations for Zg as a function 
of (j/o,2o); where yo evolves from 1 to -20(521.1,0) and zq from 20(521.1,0) to 20(521.1,521/2), equation ()5.38|l 
becomes 

(5.43) MCI^] =ZgGg+Fg, 

\dzo / 
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where 

V C^2-2o^(l - ViC^^S2i,iZ^^~^ - V2C^^S2v2ZQ'~^y^ ) ' 

and 521^1 and 82^2 ^-''^ given by equations (|5.41|l and (|5.42(l . 

If MC is invertible then Zg can be found by integrating the differential equation: 

(5.44) = [C-^M-^G,], + [C-^M-iF,], , 

where [V]2 denotes the second component of the vector V. The initial condition is that Zg{yQ,ZQ = y^) 
agrees with Theorem 12.31 with ly = ly'i. 

We will now prove Theorem 15. 51 Start by expanding the LHS of H5.34|l : 

(5.45) 

(2 - 2.g)(l - 2g)eg - {v^ - l)(i^i + 2 - '^g)s2v^egs2,.^ ~ {vi - l)(i^2 + 2 - ^g)s2u2&gs2,.2 

+ (1^1 - \ fs\^^egs2^^s2^^ + 2{v-Y - \){v2 - l)s2^iS2i/2egs2„iS2.2 + {^1 - ^fsiv2^gs2u2s^u2 = driverSg(s2i.i,S2^2). 
Changing variables to 

{S2v^,u) = (s2i/i,S2,.2S2ifr"^^^^'''"^^) 

induces 

Equation H5.45|) may be rewritten in these new variables as: 

whose LHS is identical to the differential equation H5.2n|l for eg{—S2vi)- 
Equation (|5.4()|) can be integrated to give 



(5.47) eg(-S2i.i,-S2i.2) = (I'l - 1) S2J1 / (S2,.J / (S2,.J 

Jo Jo 

. drivers,(4'.„ z^^l^r'^^''''"'^)'^^'.,^^., 

+ if,(u)4r')/(--) + x2(u)4t7^^)/(--^ 

where u = 82^2^21^1^ ^^/(''i ajj^;! i^i(u) and K2{u) are functions of u only. 
Switch the order of integration in H5.47II . 

ir. AO\ I \ I n-2 -(2-2g)/(i.i-l) /■"^"i C^"^ ( I \(2-''i)/(''i-l) / // N (2-i^i -2g)/(i/i - 1) 

(5.48) eg(-S2^i,-S2,.J = (J^l - 1) //V / / / (^g^^J (^g^^J 

JO Js" 

•drivers,(4;^,^.s"(--'^^'""'V4.,d4'., 



ifi(.)4l7^^^/^''^-^^+i^2(^.)4l7^^^^^''^-^ 



and carry out the s'^^^ integral in (|5.48|) to find 

eg (-S2,yi, -32,^2) = {Vx - 



.21]-'^^ I ^"^-'^ r (,^^J(2— 2.)/('^.-i) drivers,(4.„..1r'^'^"^'V4.2 
Jo 



_,-(2-2g)/(.-i) - (4;j(---^)/(--)drivers,(47,,...4r'''^""'V4. 

+AMH)4l:'^^^^''^-^^+i^2(^.)4l7^^^^^''^-^^. 

This proves formula H5.35|l of Theorem 15.51 
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6. Conclusions 



In this paper we have made a detailed study of the coefRcients eg{x,t) in the asymptotic expansion of 
the logarithm of the random matrix partition function (|1.7() for a single non-trivial time parameter, t = t2v 
and X = k/N near 1. These analytic coefficients are generating functions for the enumeration of g-maps. In 
particular, 

(1) We derived a hierarchy of differential equations for these generating functions. 

(2) We described a procedure for solving these differential equations. 

(3) Along the way to deriving the hierarchy we also derive an hierarchy of forced Burgers type equations 
for the auxiliary coefficients which we denote Zg(s); these are of combinatorial interest in their own 
right in that 

-g—^Zg{0) {two-legged g-maps with n 2i/-valent vertices }. 

A leg is an edge emerging from a univalent vertex; so that the leg is the only edge incident to that 
vertex. 

(4) We have calculated explicit formulae for eg{—s) for small values of The s-derivatives, of sufficiently 
large order, turn out to be rational functions of the endpoints squared, = Azq. The endpoints 
referred to here are the endpoints of the support of the associated equilibrium measure. Moreover, 
the coefficients appearing in these expressions are rational constants. 

The random matrix partition functions and their relations to graphical enumeration through diagrammatic 
expansion offer powerful tools for bringing methods of complex analysis to bear on fundamental questions in 
diverse fields such as statistical mechanics ^3 and combinatorics [261 15]. However, real progress in this regard 
has been hampered by the fact that, up till now, these connections have been based on formal procedures 
and conjectures. 

The results in this paper place many of these connections, for the large N expansion of the UE partition 
function, on a rigorous foundation. Moreover, the methods presented here, based on continuum limits of 
the Toda Lattice hierarchy, have yielded novel and effective procedures for explicitly calculating the relevant 
asymptotic generating functions. We hope that these results will help to spur renewed application of complex 
analytic methods in problems of statistical physics and combinatorics. Some results in this direction will 
appear in future work [^1^ . 

These analytical tools also point the way to novel combinatorial reults. The calculations mentioned in 
the last item of the list above, strongly motivate the 

Conjecture: The s-derivatives of sufhciently high order of the generating functions, e^, for fixed genus 
with arbitrary vertex valence number, can be expressed as rational functions of the endpoints of support of 
the equilibrium measure. Moreover, the coefficients appearing in these expressions are rational constants. 

As far as we know, such a conjecture has not appeared in the literature on combinatorics of maps. 

Another manuscript in preparation jl6| rigorously establishes a nonlocal representation for the eg(I,i). 
In the physics literature such representations are referred to as loop equations . We expect this to provide 
elegant methods that, together with the present manuscript, could enable us to prove the previous conjecture. 
Moreover, the generalizations of this loop equation together with the results of section IB. Ill can help to guide 
the characterization of eg(t) as a function of multiple times. The derivation of explicit closed form expressions 
for eg(t) will have relevance to a number of other current research programs in the statistical mechanics of 
combinatorial analysis such as random graphs, random tilings and polynuclear growth models. 
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Appendix A. 

In this appendix we show how the expansion H3.25|) for the integral 13.24|l was derived. 
Define 



(A.l) 

then IpOljl is 



Integration by parts ( u = s^'^ ^ and dv — sy's^-fP) of dj gives the equation 

(A.2) d, - ^X'^-'iX' - VA^3^_ f I^l_l)52,-2(^2 _ ^2)y^?3^^^_ 



Equation l|A.2|l produces a recursion relation for dj: 

1 1 



(A.3) 



di = — 
' 2j + l 



2 j + l 

The initial condition of this recursion is 
(A.4) do = r VA^^dA - ^A VA^^ - \o: 



■j-i- 



- + yA^^ 



The recursion relation l)A.3p . I|A.4|I is solved by 



(A.5) 



+ S\'\X) 



+sfiX)log 



X , v/A^^ 



(3 (3 



(2) (3) 

The functions S) ' and 5]^ are simple. We find the expressions: 



(A.6) 
(A.7) 



^(2) 



1\ 1 



-A" = ^A 



5", 



(3) 



^/2j 



J \2 



4J + 1 



The case of S'j'^^ is more complicated; after some work we have 



(A.8) 



-2i+l 



(2^-2^-1) 2 ^t,.(* + l)- 
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The integral ^1^, in terms of 5"]^^ S^^\ and S*]^^ is 



(A.ll) 



(A.9) / h{s)^/s^ - (3^ds 

(A.IO) 



/3 



Using the relation 13.1411 and the explicit expressions for the S'j*^'' derived above, we finally arrive at the form 
of the expansion for (|3.24|) given in (|3.25|) . 



References 

1. M. Abramowitz and I. Stogun. Handbook of Mathematical Functions, U.S. Govt. Printing Office, Washington, D.C., 1972. 

2. J. Ambjorn, L. Chekhov, C.F. Kristjansen, Yu. Makeenko, Matrix Model Calculations Beyond the Spherical Limit Nuclear 
Physics B 404 (1993), 127-172. 

3. S. Albeverio, L. Pastur, and M. Shcherbina. On the 1/n expansion for some unitary invariant ensembles of random matrices. 
Comm. Math. Phys. 224 (2001), 271-305. 

4. M. Bauer and C. Itzykson. Triangulations, Discr. Math. 156 (1996) 29-81. 

5. D. Bessis, X. Itzykson, and J.B. Zuber. Quantum Field Theory Techniques in Graphical Enumeration. Adv. Appl. Math. 1 
(1980) 109-157. 

6. A. Bloch, F. Golse, A. Uribe. Dispersionless Toda and Toeplitz operators. Duke Math. J. 117 (2003), 157-196. 

7. J. Bouttier, P. DiFrancesco and E. Guitter. Census of Planar Maps: i,From the One-Matrix Model Solution to a Combinatorial 
Proof. arXiv;cond-mat/0207682 

8. P. Bleher and A. Its. Asymptotics of the partition function of a random matrix model. Ann. Inst. Fourier (Grenoble) 55 
(2005), no. 6, 1943-2000. 

9. M. Bousquet-Melou and G. Schaeffer. The Degree Distribution in Bipartite Planar Maps: Applications to the Ising Model. 
arXiv:math.CO/0211070 

10. P. Deift, T. Kriecherbauer and K. T-R McLaughlin. New results on the equilibrium measure for logarithmic potentials in 
the presence of an external field. J. Approx. Thry., 95 (1998), 388-475. 

11. P. Deift, T. Kriecherbauer, K.T.-R. McLaughlin, S. Venakides, and X. Zhou. Uniform asymptotics for polynomials or- 
thogonal with respect to varying exponential weights and applications to universality questions in random matrix theory. 
Commun. Pure Appl. Math. 52 (11) (1999) 1335-1425. 

12. P. Deift and K. T-R McLaughlin. A Continuum Limit of the Toda Lattice. Memoirs of the AMS, 131, no. 624, January 
1998. 

13. P. Di Francesco, P. Ginsparg and J. Zinn-Justin. 2D gravity and random matrices. Physics Reports 254 (1995) 1-133. 

14. N.M. Ercolani, H. Flaschka and S. Singer. The Geometry of the Full Toda Lattice. Progress in Mathematics 115 (1993) 
181-226. 

15. N.M. Ercolani and K. D. T-R McLaughlin. Asymptotics of the partition function for random matrices via Riemann-Hilbert 
techniques, and applications to graphical enumeration. Internat. Math. Research Notices 14, 755-820 (2003). 

16. N.M. Ercolani and K. D. T-R McLaughlin. Loop Equations for the Large N Limit of Random Matrices. In Preparation. 

17. A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Higher Transcendental Functions, Vol. 1. New York: Krieger, 
30-31, 1981. 

18. H. Flaschka. Integrable Systems and Torus Actions, in Lectures on Integrable Systems, eds. O. Babelon, P. Cartier and Y. 
Kosmann-Schwarzbach. World Scientific, London, 43-101, 1994. 

19. K. Johansson. On fluctuations of eigenvalues of random hermitian matrices. Duke Mathematics Journal 91 (1998), no. 1, 
151-204. 

20. M.L. Mehta. Random Matrices, 2nd Edition, Academic Press, San Diego, CA, 1991. 

21. V. Pierce. The asymtptoic expansion of the partition function of random matrices. PhD Thesis, University of Arizona, 2004. 

22. V. Pierce. Combinatoric Results to Planar Maps. In Preparation. 

23. V. Pierce. An Algorithm for Map Enumeration. In Preparation. 

24. E. B. Saff and V. Totik. Logarithmic Potentials with External Fields. New York: Springer- Verlag, 1997. 

25. G. Szego, Orthogonal Polynomials, AMS Colloquium Publications, Vol. 23, 1939. 

40 



26. E. Wittcn. Two-Dimensional Gravity and Intersection Theory on Moduli Space. Surveys in Differential Geometry, vol 1, 
(1991), 243-310. 

27. A. Zvonkin. How to Draw a Group. Discrete Math., vol 180, (1998) 403-413. 

Dept. of Math, Univ. of Arizona, 520-621-2713, FAX: 520-626-5186 
E-mail address: ercolani@math.arizona.edu 

Dept. of Math., Univ. of Arizona 
E-mail address: mciamath.arizona.edu 

Dept. of Math., The Ohio State University 
E-mail address: vpierce@math.ohio-state.edu 



41 



